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Polarization observables in inclusive and exclusive electrodisintegration of the deuteron using a 
polarized beam and an oriented target are systematically surveyed using the standard nonrelativistic 
^> . framework of nuclear theory but with leading order relativistic contributions included. The structure 

■ functions and the asymmetries corresponding to the various nucleon polarization components are 

' studied in a variety of kinematic regions with respect to their sensitivity to realistic iViV-potential 

£vQ , models, to subnuclear degrees of freedom in terms of meson exchange currents, isobar configurations 

and to relativistic effects in different kinematical regions, serving as a benchmark for a test of present 
standard nuclear theory with effective degrees of freedom. 
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I. INTRODUCTION 



Over the past decade we have made a systematic study of inclusive and exclusive deuteron electrodisintegration 
with special emphasis on polarization observables [1-6]. The main purpose of this study was to reveal to what 
extent the use of polarized electrons, polarized targets and polarization analysis of the outgoing nucleons will allow a 
more thorough and more detailed investigation of the dynamical features of the two-nucleon system than is possible 
without the use of polarization degrees of freedom (d.o.f.). Specifically, our interest was focused on the role of the 
iViV-interaction model, of subnuclear degrees of freedom in terms of meson and isobar d.o.f. and in some cases on the 
role of relativistic effects. 

In the first paper of this investigation [1] we have considered the inclusive process, followed in [2] by the exclusive 
case, d(e,e'N)N, including beam and target polarization but without analysis of the outgoing nucleon polarization. 
In parallel, as an extension to previous work in photodisintegration [7, 8], we have formally derived in [3] all possible 
polarization structure functions, in total 648, and linear relations between them since only 324 can be linearly inde- 
pendent considering the fact that each structure function is a hermitean form of 18 independent complex f-matrix 
elements, provided parity conservation holds. Formal expressions for polarization observables, using a different rep- 
resentation scheme for the structure functions, have been given by Dmitrasinovic and GROSS [9] where also the 
question of necessary and sufficient measurements for a complete determination of all transition amplitudes has been 
discussed in detail. In [4] we have continued our own study by looking at the polarization of one or both of the final 
state nucleons in the exclusive processes d(e, e'N)N and d(e, e'NN) with various combinations of beam and target 
polarizations. Although one has 324 linearly independent observables, they are not independent in the more general 
sense of considering them as functions of the complex f-matrix elements. In view of the fact that for the 18 complex 
matrix elements one phase can be arbitrarily chosen and thus all observables are functions of 35 independent and 
real variables, it is obvious that the maximal number of independent observables is 35. Indeed, there exist quadratic 
relations among them, reducing the 324 linearly independent observables to the required number. The remaining 
question then is, which one of the many possible subsets of 35 observables constitutes an independent set. This 
question has been investigated recently in [5] for a two-body reaction of the type a + b — ► c + d, where we have 
derived a general criterion for the selection of a complete set of independent observables, which subsequently has been 
applied in [6] to the electromagnetic deuteron break-up. In the latter work we have also derived the relations between 
our structure functions and the ones of [9]. Most recently, we have obtained the general multipole expansion of the 
structure functions in [10]. 

Besides this interest in the hadron dynamics, there is a second aspect which underlines the important role of the 
deuteron from a different perspective, namely in providing an effective neutron target because of its extremely weak 
binding. In fact, quasifree reactions on the deuteron are frequently used in order to investigate properties of the 
neutron. A prominent example is the current interest in the determination of the electric form factor of the neutron 
Geu in quasifree electron scattering off the deuteron using longitudinally polarized electrons and either a vector 
polarized target or measuring the polarization of the outgoing neutron. Therefore, we also have studied in detail the 
sensitivity of polarization observables to Geu- 
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FIG. 1: Geometry of exclusive electron-deuteron scattering with polarized electrons and an oriented deuteron target. The 
relative np-momentum, denoted by p np , is characterized by angles 9 — 9 np and <j> — <f) np where the deuteron orientation axis, 
denoted by d, is specified by angles 9d and 4>d- 



The special and very fundamental role of the two-nucleon system for the investigation of the hadronic structure 
of nuclei, playing the same role as the hydrogen atom in atomic physics, is underlined first of all by the fact, that 
iViV-scattering is of crucial importance for fitting realistic AW-potential models. Secondly, the deuteron constitutes 
the simplest nucleus. It is very weakly bound and allows an exact theoretical treatment, at least in the nonrelativistic 
regime. Thus it should be clear that such an extensive survey is justified. Moreover, the electromagnetic probe allows 
a particularly clean and simple interpretation of the associated observables, because it possesses a well known but 
weak interaction so that in most cases lowest order approaches are sufficient. 

With the present work we want to give a concise and self-contained summary of this extensive study in order to 
provide the interested experimentalist and theorist a sort of handbook for the powerful tool of polarization observables 
in this fundamental process. We furthermore want to update our previous results with respect to (i) the recent high 
precision TV TV-potentials, (ii) consistent 7t- and p-like meson exchange currents, and (iii) complete and consistent 
inclusion of leading order relativistic contributions in a w/c-expansion. Thus we consider this work also as a benchmark 
for the status of present standard nuclear theory with effective degrees of freedom in terms of nucleon, meson and 
isobar d.o.f. To this end we will collect in the next section all material relevant for the formal aspects of this process 
including kinematic properties and relations, and definition of observables for the inclusive as well as for the exclusive 
reaction in terms of form factors and structure functions, respectively. In Section III we will discuss the experimental 
separation of structure functions and the question of how to select a complete set of observables. Some calculational 
details will be presented in Section IV with respect to the hadronic interaction and the electromagnetic current. The 
latter comprises contributions from one-body and meson exchange currents (MEC) and isobar configurations (IC) 
as well as leading order relativistic contributions (RC). Then we will discuss in Section V the influence of the above 
mentioned dynamical effects on various inclusive and exclusive observables for different kinematic regions of energy 
and momentum transfers. Finally, we will close with a summary and an outlook. 



II. GENERAL FORMALISM 



A. Kinematics 



At first we will consider the kinematic properties of the disintegration process e + d — ► e' + n + p. The kinematics 
is governed by the four momenta of the participating particles, i.e., the four- momenta k\ and ki of incoming and 
scattered electron, respectively, and the four- momenta pd, p p and p n of deuteron and outgoing proton and neutron, 
respectively. In view of on-shell conditions for the participating particles, four-momentum conservation, freedom in 
the choice of the incoming electron direction and in the choice of orientation of the scattering plane, one is left with 
five independent variables for this process. Different choices are possible and will be specified below. 

The scattering geometry is illustrated in Fig. 1. We distinguish three different planes which all intersect in one line 
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as defined by the direction of the three-momentum transfer q = k\ — k 2 . First, there is the scattering plane which 
is defined by the three-momenta of incoming and scattered electrons k\ and k 2 , respectively, then the reaction plane 
defined by the momentum transfer and the relative momentum p np = (p p — p n )/2 of the two outgoing nuclcons, and 
finally the orientation plane as defined again by the momentum transfer and the axis of orientation for a polarized 
dcuteron. 

The principal frames of reference are associated with the scattering plane, namely the laboratory frame, and the 
cm. frame of the final two nucleons, which is related to the former one by a boost along q. The z-axis is chosen in 
both frames along q and the j/-axis in the direction of k\ x k 2 , i.e., perpendicular to the scattering plane. Finally, 
the £-axis is defined by e x = e y x e z in order to form a right-handed system. The laboratory frame is the natural 
choice for the experimental determination of observables, whereas the cm. frame is very convenient for the theoretical 
calculation. Where necessary, we will indicate by a superscript "lab" or "cm." to which frame a given quantity refers. 
With respect to the cm. frame, we will denote throughout this paper by 9 and the spherical angles of the relative 
momentum p np = (p cm , 0, </>). Thus the spherical angles of proton and neutron momenta in this frame are 9 p m - = 8, 
(j)^- m - = (j) and 6n m ' = tt — 9, (j) c n m - = <f> + ir. The final hadronic state is furthermore characterized by the excitation 
energy E np which is related to its invariant mass W np by 

E np = W np -2M, (1) 

where M denotes the average nucleon mass. Finally, 9 d and <j) d denote the lab frame spherical angles of the deuteron 
orientation axis in case a polarized deuteron target is employed. 

The relevant quantities in the lab frame are the four-momenta of the incoming and scattered electrons fc' ab = 
(i?} ab , fc' ab ) and fcif b = (E 2 ah , ^2 ab )j respectively, the scattering angle # e ab , and the proton and neutron three-momenta 
p p ab = (p l p b , 9 l p h , 4> l p h ) and p^ ab = (pjf b , 6>Jf b , <^f b ), respectively. For the five independent variables one convenient 
choice is E\ &h , -E 2 ab , #g ab , 9 p ah , and </>p ab . All other quantities are then determined by them. An alternative choice is 
E l i h , E\ ab , #g ab , and the spherical angles of the relative np-momentum in the cm. system, 9, and <fi. Still another 
useful choice is the np- final state excitation energy E np , the three momentum transfer q c m , again the angles 9 and 
</>, all with respect to the cm. frame of the final np-state, and the lab electron scattering angle # e ab . This choice is 
particularly useful for formal investigations of the structure functions, because the latter depend solely on E np , q ™ 1 -, 
and 9, if calculated in the cm. frame. 

We will now list the various relevant kinematic quantities in the lab and cm. frames and their relation to the chosen 
independent variables. Throughout this work we will denote the square of a four-vector x v by = x$ — x 2 and 
use x = \x\. For given electron momenta kl ah , k 2 Uh and scattering angle 6> e ab , one has for the energy and momentum 
transfers in the lab frame 

^lab = £lab_£la bj ^ 
^lab = fclab_^lab j (3) 

and from this the invariant mass W np of the np-final state in terms of lab and cm. frame quantities 



sj{E^Y-{q^Y (4a) 

= y/(M d (M d + 2u 1 * b )+ql) (4b) 

= ^+£r, (4c) 

where 

E c d m - = yl Ml + (<f- -) 2 (5) 

denotes the deuteron cm. energy, and 

E l n f - c lab + M d (6) 

the lab energy of the final hadronic state. At the photon point o> lab = q c m - one finds according to (4c) and (5) the 
relation 

WL - M J 

^c.m. _ "/■ 



<*- + '->('- 2(1;;+ 2M) )' < ? » 
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where e d denotes the deuteron binding energy. 

The boost parameter 7, which governs the transformation from the lab to the cm. frame and vice versa, is given 

by 

? = (8a) 



M d 

Energy and momentum transfers in the lab and cm. frames are related to each other by 



Similarly, one has for the relative np-momentum 



„lab = 



_ _gS 

y/M% + (g c - m -) 2 cos 2 8 



p cm - = p% 



(£) 2cos2 ^ b < (ioc) 



^Wl p + {q^Y cos 2 e p b 



(8b) 



(M d ^ b + ql), (9a) 



q c.m. = i^Lg lab . (9b) 



^• m V i+ (Q 2cos2 "' (ioa) 



cosflff = =4= = cosfl, (10b) 



COS0 = ; = " P = = : COSflffi ■ (lOd) 



In the cm. frame, the nucleon energies arc (neglecting the small proton-neutron mass difference) 

E c.m. = E c.m. = E c.m. = ^np = ^ M 2 + (pC.m.)2 , (11) 

and their three-momenta are given by the relative np-momentum 

i%- m - = (p c - m ;0,<l>), Pn m - = (p c - m -,Tr-6,<b + n). (12) 
The same quantities in the lab frame may be expressed in terms of the cm. variables by 
lab E° d - m (W n p rf " cos x 

, zpc.m ~c.m. \ 2 1 / r -i \ 

(P P/ b " )2 = ( d M d ) +J^ d (( < l C - m f[ M2 + (P C - m f cos20 \ ±E d m W np q c ^ (13b) 
The Jacobian for the transformation f^™' — ► ri- ab (i € {n,p}) is given by 

afi C ' m - 1 / fllab^lab .3. o . _1 

«35r-i(^4_)(i + ^r<-«f-) ■ <») 

where the angle 0?- m - or #] ab denotes the angle of the momentum of the particle "i" in the indicated frames of reference, 
and /3 = a/7 2 — I/7. Furthermore, 

is the boost parameter that takes particle "i" from its rest system to the cm. frame. It is the same for both particles. 
Similarly, for the boost from the particle rest frame to the lab one has 

7 lab = 7 7 c.m. ( 1+ p pcm. cos gc.m.) _ 
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Note the relations 6^ m ' = 0, and 0£- m - = tt — 8 (see Fig. 1). Furthermore, the particle lab angle is given by 



s/ c.m. 

?] ah = arcsin ' 



/Qlab ~iab 



sin( 



(17) 



For the description of the polarization components of the outgoing particle, one associates with each particle "i" a 
frame of reference according to the Madison convention, for which the z-axis is taken along the particle momentum, 
i.e., in the reaction plane, the y-axis along <f x pi, i.e., perpendicular to the reaction plane, and the x-axis is then 
determined by the requirement to form a right-handed system. Often the polarization components are evaluated in 
the cm. system whereas the experimental measurement is done in the lab frame. Then it is necessary to convert these 
observables to the laboratory system. Applying nonrelativistic kinematics, the spin eigenstates in either system are 
simply related by a rotation, 9f m - — #| ab about the y-axis. However, it is well known (cf. refs. [11-13], for example) 
that for rclativistic kinematics there is a correction such that the actual angle of rotation, the Thomas- Wigner angle 
$Y i is given by 



iW 

h = arcsin 



One readily observes that for nonrelativistic boosts the Wigner angle becomes simply the 0f m - — #' ab . Since the 
rotation is about the y-axis, the y-components of the polarization of the outgoing nucleons undergo no change while 
the x- and z-components mix according to 

Pl* h (i) = R%(i)Pt m -, ^Gfen}, (19) 

where 

/ cos6f sm6f\ 
R w {i) = 1 . (20) 
\-sin0f cosef] 

Similarly, double polarization observables transform as 

^ = ^(p)^(n)P^. (21) 



B. Definition of observables 



The most general form of an observable in deuteron electrodisintegration is 



0(Q X ) = P X S 

= tr{T^ x T Pi ), (22) 



where 



d 3 a 

dk l i h dWf>dn^ K ' 

denotes the unpolarized cross section. Ox is an operator in the final two-nucleon spin space with P x as corresponding 
polarization observable. Its specific form depends on the analysis of the hadronic final state, i.e., whether or not 
polarization components of one or both outgoing nucleons are measured, and is defined below. Polarization analysis 
of the scattered electron is not considered here. T denotes the reaction matrix, and pi the density matrix for the spin 
degrees of the initial system. The trace refers to all initial state spin degrees of freedom comprising incoming electron 
and target deuteron. 

In the one-photon-exchange approximation the reaction matrix T separates into a leptonic and a hadronic part, 
and one obtains from (22) the well-known expression 

0(fl x ) = 3c(k[ ah ,k l 2 ah )tr(T^ x T Pl ), (24) 

where the hadronic part is represented by the T-matrix which is related to the current matrix element between the 
initial deuteron state and the final np-scattering state. The electron kinematics refers to the lab frame while the T- 
matrix and all quantities of the final np-state refer according to our choice to the final state cm. system. In (24) the 
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initial state density matrix p^ refers now to the spin degrees of the exchanged virtual photon and the deuteron, i.e. the 
virtual photon polarizations A (= 0, ±1) and the deuteron spin projections X d with respect to a chosen quantization 
axis, here parallel to q. Furthermore, the kinematic factor in (24) is 

c (jfeW> ,k l i h ) = A, -R— r, (25) 

with a denoting the fine structure constant and Q 2 = —q 2 . the four-momentum transfer squared (q = k\ — &2). This 
factor is related to the Mott cross section ayiott by 

/,iab ,iabx 1 tan 2 (0' ab /2) , a 2 cos 2 (#' ab /2) 

C(fc ' k2 ) = 6^ Q 2 ^Mott , where ^ t ^-^. (26) 

The explicit form of the unpolarized cross section in terms of structure functions is given below in (53). The spin 
degrees of the final state may be taken as s, the total spin of the np-final state, and m s its projection on the relative 
nrj-momentum p np in the final np-c.m. system. Another convenient choice are the helicities \ p and A„ of proton and 
neutron, respectively. 

Then the T-matrix of (24) between the initial deuteron state |A<j) and the final rip-scattering state \m\m2), both 
in non-covariant normalization, is given by 

T mimaAA <(M) = -n^2a Pnp Ec-™-E c d ™-/M d ( mi m 2 \J x (q)\\ d ) (27a) 
= e^ x+x ^t mim2XXd (e) , (27b) 

where A = 0, ±1, and the spherical angles of the relative momentum p np of the final neutron-proton state in the cm. 
system are denoted by (6,(j>) as already defined above. Here, Jo(q) denotes the Fourier component of the charge 
density operator and J±i(q) the Fourier components of the transverse current density operator. Furthermore, 
(7711,7712) stands for the spin quantum numbers of the final two-nucleon state, either in the standard (coupled) 
representation (s, m s ) of the total spin s of the outgoing nucleons and its projection m s on the relative momentum, 
or in the hclicity (uncoupled) basis (A p , A„). The transformation from one representation to the other is simply given 
by a Clebsch-Gordan coefficient 



(28) 



In Eq. (27a) noncovariant state normalization has been assumed and the hadronic cm. motion has been eliminated 
already. Thus initial and final hadronic states refer to the relative two-body motion in the hadronic rest frame. 
Eq. (27b) defines the reduced t-matrix. If parity is conserved, it obeys the symmetry relation 

ts — m s — \—\ d = (~ ) S d tsm s \\ d 1 (29) 

for the standard representation, and 

t-A p -A„-A-A d = (-) Ap+A " +A+Ad tA p A„AA d , (30) 

for the helicity representation. This relation reduces the number of independent t-matrix elements to 18, six for the 
longitudinal (A = 0) and twelve for the transverse (A = ±1) matrix elements. 

The initial state density matrix pi in (24) is a direct product of the density matrices p 1 of the virtual photon and 
p d of the deuteron 

Pi = P ^® P d . (31) 

For the evaluation of p 1 of the virtual photon, we allow the incoming electrons to be partially longitudinally polarized 
of degree h. This restriction does not mean a loss of generality because, as has been shown in [3], one obtains already 
for this case the maximal number of linearly independent observables. The virtual photon density matrix can be split 
into an unpolarized and a polarized part 



P\\> = P°\\> + Vaa- , 



(32) 
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where p° and p' can be expanded in terms of independent components p a and p' a (a <G {L, T, LT, TT}) according to 
the various combinations of longitudinal and transverse polarization 



plx = E ^\\'P a > ( 33a ) 

a=L,T,LT, TT 

Pxx' = E ( 33b ) 



a=L,T,LT, TT 



with 



^AA' 


— <5aa"5ao , 


SjLT . 
XX' ' 


= A'(5 A o + A^a'o , 


^AA' 


= 5av|A|, 


rTT . 

AA' - 


= <5a,-a'|A|, 


r) /L 
XX' 


= 0, 


A' LT 
XX' 


= |A'|J A o + |A|^A'o 


sj/T 
AA' 


= SxyX, 


r/ TT 
AA' 


= 0. 



(34) 



They obey the symmetries 

Sx'x = ^aa' = R A+A >a'-a , (35a) 
S'x'x = ^AA' = (-) 1+A+V n'-A- (35b) 
The independent components p a and p' a are given by the well-known expressions (note Q 2 = —q 2 > 0) 

PL = P°oo = P 2 Q% , PT = P°n = \Q 2 (l + , 



PLT = pgl = /5Q 2 f y ^ > ^T - P°_H = -O 2 ^ , (36) 



with 



p'lt — Poi — 3 ^yfey ^ ' Pt — Pii — 3Q 2 

n lab f\l /llab 

/3 = V' e= (^T' ?? = tan2( ^ ) ' (37) 

where /3 expresses the boost from the lab system to the frame in which the hadronic current is evaluated and q c 
denotes the momentum transfer in this frame. If, as is the case here, one calculates the observables in the final 
np-c.m. system, one has q c = q c - m - . We further note the simple relation to the often used parametrization of the 
virtual photon density matrix in terms of the quantities v a Q) of Ref. [14] (for (3 = 1) 

p£ = f^w, (38) 

where a e {L, T, LT, TT}. 

Furthermore, the deuteron density matrix p d can be expressed in terms of irreducible spin operators with 
respect to the deuteron spin space 

Ptx d ' = \Y.^) Mi ^\r [ ^ d )Pt-M, (39) 

IM 

where P ( f Q = 1, and P\_m arL d P^-m describe vector and tensor polarization components of the deuteron, respectively. 
We use throughout the notation I = \/2I + 1. The spin operators are defined by their reduced matrix elements 

(IHt^HI) = VSI for 7 = 0,1,2. (40) 

From now on we will assume that the deuteron density matrix is diagonal with respect to an orientation axis d having 
spherical angles (64, <j>d) with respect to the coordinate system associated with the scattering plane in the lab frame 
(see Fig. 1). Then one has with respect to d as quantization axis 

Pram' Prn^mm' > (41) 
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where p m denotes the probability for finding a deuteron spin projection m on the orientation axis. With respect to 
this axis one has Pf M {d) = Pf S M _ , where the orientation parameters Pf are related to the p m by 



pf = VTsiyi-y-^f 1 1 f) 



= £i,o + y^(pi-P-i)5/,i + -^(l-3p )<5/,2. (42) 

The polarization components in the chosen lab frame are obtained from the Pf by a rotation 

Pf M (z) = Pfe iM ^d I M0 (6 d ), (43) 
where d 3 mm , denotes a small rotation matrix [15]. Thus the deuteron density matrix becomes finally 

pL v - -in 1- * - E ' (i _\ , a 7 ,") pfe-^dum ■ (44) 



^ y w -a- m 

This means, the deuteron target is characterized by four parameters, namely the vector and tensor polarizations Pf 
and > respectively, and by the orientation angles #d and 0<j. Note that the deuteron density matrix undergoes no 
change in the transformation from the lab to the cm. system, since the boost to the cm. system is collinear with the 
deuteron quantization axis [16]. 

Now we turn to the definition of the operator Qx characterizing the various observablcs. One has 16 independent 
observables according to all combinations of the four operators (1 2) a) m the spin space of each of the two nuclcons. 
In detail, if no polarization analysis of the outgoing nucleons is performed, one has 

ill = ^00 = (T Q (p) <E> <r (n) , (45) 

where we have defined ao = 1-2. If the polarization component Xi of the outgoing proton or neutron, respectively, is 
measured, the corresponding operator is 

fi»o = cr»(p) ® c (ra) or O i = cr (p) ® <r, (n) . (46) 

Finally, the combined measurement of the polarization components X{ (p) and Xj (n) of both final particles is represented 

by 

Hij = <y t (p) ® Vj(n) . (47) 

Thus each observable X is represented by a pair X = (a' a) with a', a = 0, . . . , 3 and related to the operator 
ila'a — &a'(p) (X 1 <J a {n). Since the T-matrix is calculated in the np-c.m. system, the spin operators of both particles 
refer to the same reference frame with z-axis parallel to p np and y-axis along q x p np , i.e., perpendicular to the reaction 
plane. Thus the polarization components of the proton are chosen according to the Madison convention while for the 
neutron the y- and z-components of P have to be reversed in order to comply with this convention. The resulting 
observables are listed in Table I and are divided into two sets, called A and B, according to their behaviour under a 
parity transformation [8]. 



TABLE I: Notation for the cartesian components of the spin observables and their division into sets A and B. 



observable 


1 


Xp 


Vv 


Zp 


x n 










set 


A 


B 


A 


B 


B 


A 


B 






observable 


XpXn 


XpUn 


XpZn 


UpXn 


VpVn 




Zp%n 


ZpUn 




set 


A 


B 


A 


B 


A 


B 


A 


B 


A 



For real photons, the photon density matrix contains transverse components only, and thus, in order to obtain the 
corresponding observables in photodisintegration, one has to make in (24) the replacements 

PL -» 0, p LT -» 0, p' LT -> 0, , 4g s 

c(fci ab , fc 2 ab ) PT -> 1/6, h c(fcl ab , fc 2 ab ) p' T -» -PJ/6, c(fcl ab , 4 ab ) PTT -» P?/6 , 

where P^ and P^J denote the degree of linear and circular photon polarization, respectively, and P ; 7 < means linear 
polarization along the x-axis while along the y-axis for P ; 7 > 0. Furthermore, P^J > or P^J < describe right or left 
handed circular polarization, respectively. 
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Structure functions 



by 



For each observable X a set of structure functions is defined as quadratic hermitean forms of the i-matrix elements 



fL M (x) 
ti M (x) 

f%±{X) 
Pr IM {X) 

pi IM± I V \ 

JLT \ a ) 



1 + Smo 
4 



die (i S ?U% 



qoim\ 



1 + Smo 
2 

1 + 5mo 
4 

1 + Smo 
4 

1 + Smo 



3?e 
3?e 



6 * (^x /M =t (-) I+M+6x - B U X U ~ M ) 

MM _|_ ^_y+M+8x, B-jj-M-M^ 



»e (i 



x 



J 



3?e 



l 1+S ? (U° X 1IM ± (_)^+M+5 X ,B W 017-M^ 



Here <5j is defined by 



5/ = (5x,b - Sn) , with fo-.s := 



for X e A 

1 for X e B 



(49a) 
(49b) 
(49c) 
(49d) 
(49e) 
(49f) 

(50) 



distinguishing the two sets of observables A and B. In the foregoing expressions, the Ws are given as bilinear hermitean 
forms in the reduced t- matrix elements, i.e., for X = {a' a) 



X'XIM 



u a> a 



tl 



yy (m[m' 2 \a a <{p)a a (n)\m 1 m 2 }t mim2 x\ d (\ d \T M ] \\' d } . 



(51) 



Although the lA's are independent of the chosen representation for the matrix elements, their explicit form in terms 
of the i-matrix elements depends certainly on the representation for the initial and final spin states. We have already 
mentioned that two conventions are in common use, the hclicity representation with spin quantum numbers (A p , A„) 
and the standard one with (s, m s ). A third representation called hybrid basis, where the quantization axis is chosen 
perpendicular to the reaction plane, was introduced in [9]. Explicit expressions are listed in the Appendix A for the 
(s, to s ^representation. More general representations, which are obtained by arbitrary rotations of the quantization 
axes of initial and final spin states, are considered in [6]. However one should keep in mind that the observables and 
thus the structure functions are independent of the representation because they are defined as traces over the spin 
degrees of freedom (see (24)). 

Note that f°°~(X), {X) and f™+(X) vanish identically for X e A and correspondingly .C+(A), fa + ( x ) 
and f™~(X) for X e B. For this reason we often use the notation f a (X), f™(X) and f™(X) instead of f° 0± (X), 
fl^{X) and /2 0± (A), respectively. 

The structure functions f$ IM ^ ± \x) (a € {L, T, LT, TT}) (primed and unprimed structure functions f' a IM ^ (X) 
and fa Mt " ± \x) are here referred to collectively as fjp™^ ± \x)) contain the complete information on the dynamical 
properties of the n/j-system available in deuteron electrodisintegration. They are functions of the np-angle 6, the 
relative np-energy E np , and the three-momentum transfer squared (g c m ) 2 , all in the cm. system. 

In terms of these structure functions a general observable in d(e, e'N)N and d(e, e'np) is given by 

2 I 

0(CI X ) = c(fcl ab , fc' ab ) E P /"E { {pLfl M (X) + PT ti M (X) + p L Tfl% + {X) COS0 



7=0 



M=0 



+PTT.fTT + (X) COS 20) CO S (M0 - Sf '-) 



- ( PltJlt ~ (X) sin 4> + pttJtt " ( X ) sin 20 ) sin(M0 - Sf - ) 



+h 



(p T fT IM (X) 



+ PLTfLT~(X) cos A sin(M0 - If J) 



+p' LT f'lT + (X) sin0cos(M0 - If |)] }d M0 (6 d ) , 



(52) 
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where we have introduced cf> = <fi — </>d- ln particular, one obtains for X = 1 and Pf = 5i$ the unpolarized cross 
section as 

S = c(fci ab , fc' ab ) ( PL f L + pt/t + PltJlt cos <f> + ptt/tt cos 20) , (53) 

using as a shorthand f a = /q 0+ (1). One should remember that the nucleon angles and polarization components refer 
to the cm. frame. The transformation to the lab frame is described in the previous subsection. 

TABLE II: Listing of (IM)- and (/Af±)-values of nonvanishing structure functions f a (X) and f a (X), respectively. 



/04 fr M (X) f?r IM (X) f I L ^ ± (X),f^ ± (X),f'^ I± (X) 



x e A 


X G B 


X G A 


X G B 


X £ A 






X € B 


00 


10 


10 


00 


00+ 


10- 


10+ 


00- 


11 


11 


11 


11 


11+ 


11- 


11+ 


11- 


20 


21 


21 


20 


20+ 


21- 


21+ 


20- 


21 


22 


22 


21 


21+ 


22- 


22+ 


21- 


22 






22 


22+ 






22- 



TABLE III: Number of nonvanishing structure functions fL^™ ' (X) and f a '^ IM± (X) for 


an observable X e A 


or X G B. 


set L T T' LT LT' 


TT 


total 


A 5 5 4 9 9 
B 4 4 5 9 9 


9 
9 


41 
40 



The possible (/, M)-values are listed in Table II and the total number of structure functions for each observable 
X and each a are listed in Table III. As mentioned in the introduction, one finds altogether a total number of 648 
observables, each of which is a hermitean form of the ^-matrix elements. However, since the ^-matrix has only n = 18 
independent complex amplitudes, only n 2 = 324 linearly independent hermitean forms can exist. Indeed, one finds n 2 
linear relations between the observables which are presented in the next section. The remaining structure functions 
are linearly independent so that indeed the maximal information can be obtained by using longitudinally polarized 
electrons alone. Transverse polarization is not necessary. 

On the other hand, since each reaction matrix element is in general a complex number, but one overall phase is 
undetermined, a set of 2n — 1 properly chosen observables should suffice to determine completely all matrix elements. 
This seeming contradiction is resolved by the observation, that the linearly independent observables are not completely 
independent of each other in a more general sense. In fact, any bilinear form t*,tj can be expressed as a linear form 
in the observables (see [5] and also Sec. II D), and for these bilinear forms one can find exactly (n — l) 2 quadratic 
relations (see Appendix B), thus reducing the total number of independent observables just to the required number 
2n — 1. Consequently, one can determine all matrix elements from 2n — 1 properly chosen observables. However, one 
should keep in mind that the solution is in general not unique but contains discrete ambiguities. This is discussed in 
Sect. III. 

To close this section, we will give for the transverse structure functions (a G {T, T', TT}) the correspondence to 
the observables in photodisintegration derived in [8]. The formal definition of observables is completely analogous 
except for the fact that in photodisintegration only transverse current components contribute. Taking into account 
the slightly different definition of the T-matrix (compare T of (27a) with the definition of T 7 in [17, 18]), i.e. 

IW nC-m. 

T mim2XXd (0,4>) = V % T ^m 2 xx d (0,<t>), (54) 
one has the following relations at the photon point with respect to the general form of an observable in photodisint- 
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TABLE IV: Definition of X'(X) and p(X). 



X 


1 




•£pZn 


y P 


Xp 




Zp 




X'(X) 
P(X) 


VpVn 
1 


Zp^n 

-l 


Zp^n 

l 


l 


^pVn 

-l 


VpZn 

-l 


^pDn 

l 


Vp^n 

l 


set 


A 


A 


A 


A 


B 


B 


B 


B 



^ration as given in (14) and (15) of 



^ = (55.) 
f T IM {X) = - 6W Z qC ' m ' Px IM , (55b) 



M d ~ x dVL 



Hp 



D. Linear relations between structure functions 

As is shown in detail in [3], the derivation of linear relations among observables is based on the inversion of (51) 
expressing any bilinear form t*, , y ,t s1Jls \ m as a linear superposition of observables. This inversion can be done 

analytically (see [3]). In general one obtains two types of relations among the structure functions of an observable X 
and those of another observable X'(X), uniquely related to X. Explicitly, one finds as the first type of equations 

g™(X) = IpiX^iX') V2f\X>) VTifHX')) , (56a) 
g n (X) = p(X)g ll (X'), (56b) 

g 2 \X) = \p{X)[-^f\X') + 2g 2 \X')-^\g 22 {X')), (56c) 
g 21 (X) = p(X)g 21 (X'), (56d) 
9 22 (X) = -^=p(X)(g m (X>) + -^9 2 °(X')) . (56e) 

The second type of equations reads 

h w (X) = --Lp(X)h 22 (X'), (57a) 

h n (X) = p(X)h 21 (X'), (57b) 
h 21 (X) = - P (X)h 11 (X'), (57c) 
h 22 (X) = V2p(X)h w (X , )) 7 (57d) 

where X'(X) and p(X) are listed in Table IV for 8 observables. For the remaining other 8 observables, not listed in 
Table IV, one obtains X'(X) from Table IV with the help of the relation 

X'(X'(X))=X, (58) 

and p(X) from the relation 

p(X'(X)) = (-)'*■* p(X). (59) 

Table V shows which of these two types of relations holds for a specific structure function, depending on whether X 
belongs to an observable of set A or B. Which structure functions are related to each other in these relations, i.e., 
g IM (X) to g I,M '(X') and h IM (X) to h I,M '(X'), is also listed in Table V. 

At the end of this section, we will give two examples of how to find the proper relation for a given structure function 
of an observable X. As first example we choose the y-component of the neutron polarization, i.e. X = y n belonging 
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TABLE V: Listing of structure functions g IM (X) and g 1 M (X 1 ) for observables X, X' which fulfill the relations (56), and of 

structure functions h IM (X) and h 1 M (X') for observables which fulfill the relations (57). The associated observable X'(X) is 
either listed in Table IV or can be obtained using (58) . 



X € A [B] 


g[h] IM (X) 
g[hf M '(X') 


flM 

Jl 

fi'M 1 

Jl 


fIM+ 
J LT 
rl' M' + 

Jlt 


rl IM+ 

Jlt 

fl i' m' + 
Jlt 


flM 

Jt 

fl'M' + 

Jtt 


fIM+ 

Jtt 

fi'M' 

— Jt 


X £B[A] 


g[h]' M (X) 
g[hf M '(X') 


flM- 

Jlt 

fll'M'- 
■1 LT 


rl IM- 

Jlt 

_fl'M'- 
J LT 


fl IM 

Jt 

fi'M' — 
JTT 


flM- 

Jtt 
ii I'm' 
~Jt 




TABLE VI: Listing of the matrix Uf s ^. 


1 


S 


A* = 1 


2 


3 


4 


i-i 


1 


COS €j 





— sin tj 





3 








1 








3 + 1 


1 


sinej 





COS tj 





j 


1 











1 



to set A. According to Table IV its structure functions are related to the ones of the y-component of the proton 
polarization. With the help of (58) and (59) one finds X'(y n ) = y p and p(y n ) = 1 and in view of Table V the relations 
(56) apply, e.g. 

fTiVn) = \ (f° L °(y P ) V2fl°(y p ) V3fl 2 (y P )) , (60) 

f T °(yn) = -\(f!&(Vp) V2f T ° T (y p ) V3ti 2 T (y p )) . (61) 

For the second example we choose X — z p y n , belonging to set B, and the structure function /^°(z p 7/„). From 
Table IV with the help of (58) and (59) one finds X'(z p y n ) — x p and p(z p y n ) = 1. Furthermore, according to Table V 
the relation (57) applies, resulting, for example, in 

fl (z p y n ) = -^=fl 2 (x p ). (62) 



E. Multipole decomposition 

A convenient parametrization of the angular dependence of observables and structure functions is provided by an 
expansion in terms of the small rotation matrices d 3 m , m [7, 19-21]. Explicit expressions for deuteron electrodisin- 
tegration have been derived recently in [10]. They facilitate the analysis of the contributions of the various charge, 
electric, and magnetic transition multipole moments to the different structure functions. This expansion is based 
on the multipole expansion of the ^-matrix. We take the outgoing np-state in the form of the Blatt-Biedenharn 
convention [22] 

\psm s )^ = Yl Hl0sm a \jm a ) e~< Ifi mjm ,[R) \m) , (63) 

Hjrrij I 

where D 3 m m (i?) denote the rotation matrices in the convention of Rose [15] and fi = 1, ... ,4 numbers the four 
possible partial waves for a given total angular momentum j > 0. For j = one has only two partial waves. The 
phase shifts are denoted by 8^, and the matrix [/j? is determined by the mixing parameters Cj as listed in Table VI. 

Furthermore, R = (0, —9, —(f)) rotates the chosen quantization axis into the direction of the relative np-momentum 
p. The partial waves 



lw'"»j> = U i's'n W s ')i m j) 

l's' 



(64) 
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TABLE VII: Listing of the sets kx determining the summation values k in the multipole decomposition (73) of a structure 
function for an observable X — (a' a). 



a' 
a 


3 3 
3 3 


1 2 3 1 3 2 
1 2 1 3 2 3 


12 1 2 
112 2 




{0} 


{-1. 1} 


{-2, 0, 2} 



are solutions of a system of coupled equations of iVTV-scattering. In this convention, the i-matrix reads 

i 



t S m s xxA0) = (-)Vi + <5ao -Al^LX\jm j )(l0am a \jm a )O LX (jijla)4 njmt (e), (65) 

Lljrrijfi 

with 

O lx {mIs) = Vtoe^U^Nfoj) , (66) 

and 

JVa («) = 5 W (E L ( M ) + \M L {^)) + 6 xa C L (w) , (67) 

where E L (/ij), M L (fij) and C L (/zj) denote the reduced electric, magnetic and charge multipole matrix elements, 
respectively, between the deuteron state and a final state partial wave \fj,j) in the Blatt-Biedenharn parametrization. 
Parity conservation implies the selection rules 

{C/E) l {m) = for (-) L +i+» = -1, (68a) 

M L ([ij) = for (-) L +i+» = l, (68b) 

which leads to the relation 

O l ~ x {mIs) = (-) L+l O LX (»jls) . (69) 

As is shown in detail in [10], one obtains for a structure function the general multipole expansion according to the 
expressions in (49) from the one oiU xx J M in (51) which reads 

jjX'XIM _ 1 ,X'XIM,Kk ,K ( n\ / 7n \ 

where d 3 ,(0) denotes the small d-function of the rotation matrices. 

tit ifi * ' 

The sets Kx of the possible k- values are listed in Table VII and the coefficients are given by 

U x : x J M ' KK = Wc-'.-, C x ' XIMK {L'j'Lj)V«l{n'j'M) N^^'f) Nfoj) , (71) 



where we have defined 



^(a'.a) = S °",k + Sa.k , (72a) 



and 



C x ' XIM > K (L'j'Lj) = (-) X '+ L 2 (1 + 5 yo )(l + 5 X0 ).fjlk 2 

».*»'/w) = h)C-> £ v$i<jwmvf,„vL,- <72i:) 
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TABLE VIII: Listing of the values of (3(a) in the multipole decomposition (73). 



a 


L/T 


LT 


TT 


(3(a) 





1 


2 



with 



E« 2 it ' ' 



1/ V —K 




L^f -k 



.S /\ /■."'• / /<■' / /' 




S K K' 
s j I 
s' f V 



(72d) 



The definition of s™ is given in (A3) of Appendix A. Furthermore, in (71) we have incorporated the phase shift for 

convenience into the quantity N^(/j,j) = e lS " N^(nj). 
Then the general multipole decomposition reads 

(73) 



K, k£kx 

where 13(a) is listed in Table VIII and the coefficients f£ IM ^' Kk (X) are obtained via (49) from the foregoing 
multipole expansion. In detail one has for the longitudinal and transverse structure functions of an observable X 



f IM, Kk/ v \ 

Jl \ a ) - 



]T C[ M,K (L'j'Lj) T>K, K a (n' j' nj) ^te(i^ +5 (^'-"i C L '* (fi'f) C L (iij)} 



L'n'j'Liij 



Jt \ a ) 



pi IM, Knz v \ 
Jt l A ) 



E C T M ' K (L'j'Lj) V«« a {pl3'M) Me(f> + V,«)< V/) Nfojj) , 
E C^ K {L'j'Lj)V^j'^rn{f'^ 



L'n'j'Lfij 



L'n'j'Lpj 

and for the interference ones, distinguishing observables of type A 



( IM±.Kki v \ 

Jtt l A ) 

£ IM±,Kk/ 



I LT 



E C™ ± ' K {L'j'Lj)l%,Wj'M)te(f? + ^ 



L'n'j'Lftj 



and observables of type B 

c IM±,Kk/ 



TT 



\X) 



rIM±. Kk/ v \ 

Jlt \ A ) 



f/IM±,K K 

Jlt 



(X) 



E ^^(Wj)^/^ 

L'n'j'Lfij 

E Cl^ K {L>j>Lj)V«r a {»>3'»me^ 



L'n'j'Lfij 



(74a) 
(74b) 
(74c) 

(75a) 
(75b) 
(75c) 

(76a) 
(76b) 
(76c) 



L'n'j'Lpj 
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Here the coefficients C a are denned by 

(77a) 
(77b) 
(77c) 
(77d) 

More detailed expressions for the coefficients of the structure functions of the differential cross section are listed in 
the Appendix C. Explicit results for the coefficients C and V for a maximal multipolarity L max — 3 may be found 
in [10]. 



Cl M > K (L>j'Lj) = 


87T 


C 00IM ' K (L'j'Lj), 




1 + Smo 




C^ K {L'j'Lj) = 


16tt 


C^M,K {L , fLjh 




1 + S M o 




C[^ K (L'fLj) = 


16 7T 


(c 01IM ' K (L'j'Lj) ± (-)i+M C au 


- M > K {L>i>Lj)) , 


1 + Smo 


C^ K {L'j'Lj) = 


8tt 


(c- 11IM ' K (L'j'Lj)±(-) I+M C- 


ni-M,K {L , fLj ^ 


1 + $M0 



F. Inclusive process and form factors 

The inclusive cross section is obtained by integration over the solid angle £1%™' — (6, <P) yielding 

dk] l a da u h = 6c ( fc 'i ab ' fc 2 ab ) {plF l + PT F T - P^ltF 1 ^ 1 sm<t> d d{ Q (e d ) 

+ P?[(plFI° + PT F^)d 2 m (9 d ) -pltF^ 1 cos fad 2 w (9 d ) + p TT F 2 T 2 cos2^ d d 2 2Q (9 d )] 
+ hPi[-p' T F™dl {6 d ) + p'ltF'It 1 cos^dio(^)] - hP^p' LT F 12 ^ sin<Mio(^)} 
= a(h,P?,P*). (78) 

It is governed by a set of inclusive form factors F^) 1 M (M > 0) as given by 

F^- M = (-Y+ M (1 + S M0 ) J I d(cos0)(/M' M + - /W JM ") , (79) 

for a e {L, T, LT, TT}. This equation corresponds to Eqs. (13) and (14) of [1] except for the fact that the primed 
form factors F^ 10 and F'jjf 1 differ in sign from the ones given in [1] due to a redefinition of the primed structure 
functions incorporating a phase factor (— ) 7 (see the remark in [4] before Eq. (9)). Altogether, the inclusive cross 



section depends on ten form factors: F L , F T , F^ 1 , F 20 , F 20 , F£^\ F 2 ^ 2 , F^ 10 , F'jjf 1 , and F 12 ? 1 , of which F^ 1 
and F'jlf 1 vanish below pion threshold due to time reversal invariancc. 

The multipole decomposition of the inclusive form factors are given by the (K ~ 0)-cocfficicnts of the multipolc 
expansion of the structure functions of the differential cross section (X = 1 = (0, 0)), as listed in (73), i.e. 

F W-M = { _y+M (1 + Smo) I (/ ( /M+, o _ ; w/m-, 0) (80) 

Explicit expressions are listed in Appendix C. 

At the photon point one can relate the purely transverse form factors to the total photoabsorption cross section a to t 
of deuteron photodisintegration for unpolarized photons and deuterons and to the corresponding beam and target 
asymmetries of the total cross section as defined in [18]. Taking into account the relations (55) one obtains respectively 

M d F F T° c F T W J F TT (RU 

= W n p q cm - ^ ' T20 = 7V Tl0 = ^T' T22 = ^~' (81) 

where W np and q c - m - denote the invariant mass of the np system and the photon cm. momentum, respectively. 



III. SEPARATION OF STRUCTURE FUNCTIONS AND COMPLETE SETS 
A. Experimental separation of structure functions 



The experimental separation of structure functions has been discussed in detail in [2, 4]. It is based on the general 
definition of asymmetries of a polarization observable with respect to the beam and target polarization parameters h, 
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Pf, and P$ , respectively. To this end one writes a general polarization observable (D(Qx) as given in (24) and (52) 
in the form 



Px = A (X) + Pt A v d {X) + P 2 d A T d {X) + h A e (X) + P* A v ed {X) + P 2 d A T ed {X) 



(82) 



defining implicitly the various asymmetries, i.e. A e (X) for beam polarization, A d (X) and A d (X) for vector and 
tensor target polarization, respectively, and A^ d (X) and A^ d (X) for the corresponding beam-target asymmetries. 
Their explicit form can be read from (52) 



A Q (X) 



A v d {X) 



c(fc^ b .fci, ab ) 
So 



c(fc 1 . fc 2 ) 
So 



AJ(X) = 



c(k 1 . k 2 ) 



(p L fT{X) + p T f^(X) + p LT fl° T + (X) cos0 + p TT f T °+(X) cob2(^Sx,a 
+ (pltIlt (X) sin + pttItt ( x ) sin 2( t>) Sx,b] , 

E [(pLfl M (X) +PT.f 1 T M {X) + PLT .fl A 4 + {X)cO^ 
M=0 

+PttJtt + {X)cos2c^ cos{M4>-5 x ,a\) 

-(pLTfl^-(X) sin0 + pTTf^-(X) sin 20) sin(M0 - <5x,a|)] <W^) , 
£ [{pLfl M {X)+ PT f T M {X)+ p LT f L A 4 + {X)co^ 



(83a) 



(83b) 



M=0 



+PTT.fTT + (X) cos 2^) cos(M0 - S x ,b\) 

PltIIt' (X) sin + PttItt' (X) sin 20 ) sin(M0 - <5 X , B 



Mi} 



A e (X) 
A v ed {X) 



So 

ab ,fc 

So 



- (/^°°P0 + P'LTf' L T-{X) COS CP)S X ,B + p' L TfLT + (X) ^5 X ,A 
1 

£ [(/t/t 1M (*) + Plt/lt M "(^) cos^) sin(M0 - ^, A |) 



M=0 



AUX) = 



_ c(fc£ 



So 



+PLT.fLT + (X) sin 0cos(M0 - 5x,a|)] 4f (^) , 
2 

E [(pt/t M P0 + p'LTfLT M -{X) cos^) sin(M0 - ^, B |) 



M=0 



Vlt/Lt (*) sin0cos(M0 - «5x,b-) 



d 1 

«M0 



(83c) 
(83d) 

(83e) 

(83f) 



where the unpolarized differential cross section So is defined in (53). For simplicity, we will also call Aq(X) an 
asymmetry although it is not one in the strict sense. The nonvanishing structure functions contributing to an 
asymmetry of a given observable are listed in Table II. For the differential cross section (X = 1 = (00)) we remind 
the reader that one has with respect to the notation in [2, 40] 



4>(1) = 1, A v d ' T {l) = A v d ' T , A e (l) = A e , A v J d T {\) = A v f . 



(84) 



For the simplest case, namely in the absence of beam and target polarization, the four structure functions f a (X) 
can be separated choosing first different 0-anglcs, yielding Jlt, Itt and a linear superposition of Jl and fx and 
subsequently a Rosenbluth separation for disentangling Jl and fx- In the general case, by a proper variation 
of the longitudinal electron polarization h and the deuteron vector and tensor polarization parameters Pf and P d , 
respectively, one can first separate the various beam, target and beam-target asymmetries as listed in (83). These 
asymmetries are functions of the deuteron orientation angles 9 d and <j> d , viz. <p — 4> — (f) d , and the azimuthal or out- 
of-planc angle <p. One can now utilize these variables for the further separation of the different structure functions. 

This is achieved by observing that the general functional form of an asymmetry is 



A 1 ^, 4>, 6 d ) = E «/m(0, 0)dk o (0«O - i 1 = °> i. 2 ) - 



(85) 



M=0 
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where 

aiM(4>, 4>) = c IM {4>) cosMcf) + s IM {4>) sin M</> , (86) 
and the (/(-dependent functions cjm(<I>) and Sim{4>) have either the form 

a + ai cos (f) + a 2 cos 2<f> (87a) 

or 

6isin0 + fe 2 sin2(?!> . (87b) 

For a given / the M-components a/M(</>, </>) °f the asymmetry A 1 ((j),(j),0d) can be separated by a proper choice of 
9d exploiting the properties of the small (ij^Q-functions. For 1=1 (vector asymmetries), taking 9 a = or it/2, i.e. 
^mo(O) = ^mo or ^0(77/2) = M/y/2, yields a\o or an, respectively, and for the tensor asymmetries (/ = 2) one may 
first choose 9d = yielding with d 2 MO (0) = 8 mo directly a 2 o- The latter being determined, then setting 9d = 7r/4 and 
7r/2, one can obtain the remaining two terms a 2 i and a 22 . For the separation of a 2 i and oj 22 one can also choose 
Sd = Q'd — arcos(l/v / 3) together with <fi and (f> + ir. Then the sum and difference of the corresponding asymmetries 
result in a 2 i and a 22 , respectively. 

In the next step, in order to separate the two contributions cim and s/m in (86), one can take first = giving 
c/m and then <p = ir/2M for M ^ which yields directly s/m- The remaining separation of the coefficients a n or b n 
in (87) is then achieved by appropriate choices of (f>. In a few cases the constant term a in (87a) will contain two 
structure functions in the combination pLf[ M (X) + pTfx M {X). In this case one needs a Rosenbluth separation in 
addition. 

A different task than the complete separation of all structure functions is to find an optimal way for separating a 
specific structure function fa M ^\x) (I > 0) or f'a IM<,± \x). In other words, the question is: what is the minimal 
number of measurements necessary for the separation of a specific structure function? This has been discussed in [4] 
and is described in detail in Appendix D. 

B. Complete sets of observables 

We have already mentioned the fact that in deuteron electrodisintegration the total number of independent complex 
i-matrix elements is 18, while for photodisintegration the number is 12. Since one phase remains arbitrary this means 
that this process is determined by 35 independent observables, whereas in the corresponding photoreaction one needs 
23. The question is how to choose from the much larger set of 324 (or 144 in photodisintegration) linearly independent 
observables an appropriate set of 35 (or 23). In [5] we had derived a general criterion which allows one to decide 
uniquely whether for a reaction with n independent ^-matrix elements a set of 2n — 1 observables, taken from the set 
of n 2 linearly independent observables, constitutes a complete set. Subsequently this criterion has been applied to 
deuteron electro- and photodisintegration in [6]. A brief review of the main results of [5] and [6] is appropriate. 

1. General criterion [5] 

Any observable in a reaction with n independent, complex matrix elements can be represented by a n x n hermitean 
form f a in the complex n-dimcnsional variable z = (z\, . . . , z n ) 

H*) = !XV'V,-, (88) 
j'j 

where hermiticity requires 

(F? j )*=F? j ,, (89) 

and z comprises all independent reaction matrix elements labeled by j. 

For the application of our criterion, derived in [5], one first has to rewrite the hermitean form in (88) into a real 
quadratic form by introducing 



z = x + iy , 
F a = A a +iB a , 



(90) 
(91) 
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where A a and B a are real matrices, and A a is symmetric whereas B a is antisymmetric. Considering further the fact 
that one overall phase is arbitrary, one may choose yj = for an arbitrary index j and then one finds for the given 
observable 



3 X 3 



(92) 



3 3 



where the tilde over a summation index indicates that the index j has to be left out. Introducing now a (m = 2n— 1)- 
dimcnsional real vector u by 



u = Oi, ...,x n ,yi,.. . ,y 3o -i,y 3o+1 , ...,y n ), 
one can represent the nx n hermitean form by a to x to real quadratic form 



i=i 



where the to x m-matrix F a is given by 



A a (B a ) T 
B a A a 



(93) 



(94) 



(95) 



Here B a is obtained from B a by canceling the jo-th row, and A a from A a by canceling the jo-th row and column. 
Thus B a is a (n — 1) x n-matrix and A a a (n — 1) x (n — l)-matrix. 

Now, for checking the completeness of a chosen set of 2n— 1 observables one has to construct the m x m corresponding 
matrices F a , and then one builds from their columns for all possible sets {k\, . . . , k m ; k a <G {!,..., to}} the matrices 



/ F 1 



W(k u ...,k m ) 



ifei 



F, 



(96) 



rnk n 



Note that the k a need not be different. If at least one of the determinants of W(ki, . . . , k m ) is nonvanishing then one 
has a complete set. 



2. Complete sets for photo- and electrodisintegration [6] 

In order to apply our criterion, one has to construct the matrices F which represent the structure functions as 
hermitean forms in the reaction matrix elements as 

f£ IM± (X) =t^F^ IM± ' a t, (97) 

where t is a vector comprising all reduced i-matrix elements in a certain labeling. It is convenient to arrange the 
labeling of the i-matrix elements in such a way that the longitudinal ones belong to j = 1 , . . . , 6 and the transverse 
ones to j = 7, . . . , 18. Thus the general structure of these matrices then is 

~, > / A(/)IM±.a r<(')IM±,ct\ 

I (CM IM± ' a )t B(')IM±,a J ' v ' 

where A^ IM± ' a is a (6 x 6)-matrix, C^ IM± > a a (6 x 12)-matrix, and B^ IM± ' a a (12 x 12)-matrix. In particular 
one has 

„ / aim±, l n \ 

piM±,L = { A Q °\, (99a) 

j?{l)IM±,T/TT _ I \ , , 

I B^ IM± - T / TT j ' 
~ ( n n(') im±, lt \ 

p(,)IM±,LT _ I U t w fg9 -> 

I ( C (/)JM±,LT)t i ' 1 ^ 
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The explicit forms of these matrices are obtained from the matrix representation of 

Ux XIM = J2 t}C I ^ yx {X)t J . (100) 
j'j 

Comparison with (51) gives for the matrix elements 

C I f f x ' x (X) = (m / 1 m 2 |<7 a ,(p)a a (n)|mim 2 )(A d |r^|A / d ) , (101) 

where the labeling is to be understood as j'W = (m^, A^). Detailed expressions of the C IMX x (X) , s for 

several representations are easily obtained from the expressions listed in Appendix E (see also [6]). 

The structure of these matrices is such that the longitudinal (L) and the transverse (T, TT) observables are decou- 
pled filling separated 6x6- and 12 x 12-submatrices, respectively, whereas the LT-type observables are represented 
by 18 x 18-matrices. These features offer various kinds of strategies for selecting complete sets. 

(i) One may independently select complete sets of observables for the longitudinal and transverse cases, i.e., a set 
of 11 longitudinal and 23 transverse structure functions for a check of completeness. With respect to the latter, 
one has in view of the linear relations between the T- and the TT+-type and between the T'- and TT^-type 
observables, different choices, taking either T- and T'-type or TT ± -type observables or even mixing different 
types of observables. The missing relative phase between the longitudinal and transverse t-matrix elements can 
then be provided by any one of the TT-observables. The advantage of this approach is that in this way one 
automatically obtains complete sets of observables for the case of photodisintegration as well, namely from the 
transverse ones. 

(ii) Again one may start with a selection of 11 longitudinal structure functions. But then instead of choosing trans- 
verse observables, one may directly choose 24 linearly independent LT-type observables which then constitute 
a simple system of linear equations for the missing transverse matrix elements, because the longitudinal ones 
are then known from the first step. 

(iii) Complementary to case (ii) one may start with a selection of 23 transverse observables taking one of the 
alternatives listed in (i). Then a proper set of 12 LT-type observables provides a set of linear equations from 
which the missing longitudinal Lmatrix elements can be obtained. 

(iv) Another alternative would be a selection of 35 structure functions of LT-type. However, in this case the 
completeness check would be much more involved due to the considerably higher dimension of the determinants 
to be checked. 

Which of these strategies is most advantageous will depend on the experimental conditions. Often L- and T-type 
structure functions arc easier to determine in an experiment although the required Rosenbluth separation introduces 
some unwanted complication. In view of the fact that the strategies (i) through (iii) require the determination of 
either L- or T-type observables or both, we have considered in [5, 6] exclusively the question of complete sets for 
longitudinal and transverse structure functions. 

The choice of complete sets of longitudinal structure functions, containing eleven structure functions, from a set of 
linearly independent ones has been discussed in detail in [5]. It turned out, that there is only a very weak restriction on 
the choice of possible complete longitudinal sets. In fact one may select from the chosen set of 36 linearly independent 
observables any subset of eleven structure functions, which docs not contain more than eight of the type X 10 and X 22 . 
This has been discussed explicitly for the linearly independent set X <G {1 /M , xx IM , xz IM , y[ M , x[ M , x^ M , z[ M , z| M }, 
and possible complete sets are listed in Tables 3 and 4 of [5]. The case of the transverse observables for which a 
complete set contains 23 structure functions, has been discussed in [6]. Again with respect to the general question of 
a choice of a complete set of structure functions, we found the general statement, that one may pick from the chosen 
set of 144 linearly independent ones any subset of 23 structure functions with the only restriction, that not more than 
16 should be of the type I« 10 and I« 22 . 

In Ref. [5] we furthermore simulated an experimental study for the determination of the longitudinal t-matrix 
elements in the helicity basis from a given set of "measured" observables whose numerical values were taken from a 
calculation. Various complete sets were selected and the arising system of 11 nonlinear equations for the t-matrix 
elements was solved numerically. Since the solutions were not unique we had to calculate additional observables, called 
"check observables", taking as input the obtained solutions for the Lmatrix elements, and compared them to their 
"measured" values. For the arbitrarily chosen kinematics (internal excitation energy E np = 100 McV, momentum 
transfer (q c m ') = 5 fm~ 2 , various np-anglcs 6) we found that one of the considered complete sets was particularly 
suitable (first set of Table 6 in Ref. [5]). In this case only one additional check observable (/i, ^)) was sufficient 
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to determine the correct solution. In [6] this simulation was extended to a somewhat more realistic experimental 
situation using the same kinematics again and taking the same specific set but allowing for errors in the measured 
observables. The simulation of an experimental situation showed that one can get quite reliable results for the i-matrix 
elements even if experimental errors are taken into account. The results can be greatly improved if additional check 
observables are considered. For the case studied in [6], it was sufficient to consider two such observables. If on the 
other hand one uses no check observables at all, one gets rather unreliable results since other types of solutions of 
the nonlinear system of equations are mixed in. In fact, performing such a simulation without any check observable 
leads to large errors in the resulting t-matrix elements (average error more than 100 %) and also to strong average 
deviations of the mean values from the true values of the i-matrix elements (about 50 %). 

A similar study has been performed in [6] for the determination of the transverse i-matrix elements from observables, 
but without introducing experimental errors. Although the transverse case is much more complicated than the 
longitudinal one due to the higher dimensionality (12 instead of 6 complex i-matrix elements), it was found that in 
principle the method works also for the transverse case. 



C. Analytic expressions of the t-matrix elements in terms of observables 

One can also derive an analytic solution of the reaction matrix elements in terms of observables, because one can 
express all bilinear forms t*,tj as linear forms in the structure functions fa^ IM± (X), i.e., 

t*, tj =T j , j [fW IM± (X)} = J2 (Tf J Msi9 fL Msi9 (X)+T^ Msi9 f a IM ^(X)), (102) 

alM.sig— ± 

where the explicitly appearing structure functions f£ IM± (X) constitute a complete, linearly independent set. Here 
the square bracket of Tyj [/] indicates the functional dependence on the structure functions. This relation has also 
been used in the derivation of the linear relations between observables (see Sect. II D). Explicit expressions for the 
coefficients j^') aIMs *9 are derived in various representations in [6]. To give an example, we list here two cases, a 
diagonal and an interference term, of the longitudinal i-matrix elements in a rotated helicity basis, where the final 
helicity states are rotated into the y-axis, 

*|i 00 *iioo = l(fl° -tiff + /E°(yi) - V2/!°(yi)) , (io3a) 

f_i-ioo*H«> = l(-fl\zz) + V2fl (zz) + i(f™(xz)-V2fl (xz))). (103b) 

For a complete listing and further details we refer to the Appendix E of [6]. 

The linear relations in (102) can be exploited in various ways. One possibility is to choose a specific matrix element, 
say tj , as real and positive. Then all other matrix elements tj with j ^ j arc uniquely determined relative to tj 
and are given as linear forms of appropriate structure functions [8] 

Finally, for the determination of the missing matrix element tj one has to choose only one additional structure 
function, say 



/o = E4^i- ( 105 ) 



yielding 



t j0 = 4=. iy.Tr^ IM± {Xm l0 T 0Q] [f^ IM± (X)\ . (106) 



Thus (104) in conjunction with (106) constitutes a nonlinear functional in the structure functions fa^ IM± (X). How- 
ever, proceeding in this way, one needs in general a much larger number of observables for the complete determination 
of the t-matrix than the required minimal number of 2n — 1 of a complete set of a n-dimensional i-matrix. 

Another strategy, which leads in general to a smaller number of necessary observables, has been developed in [6]. It 
is based on an analysis of all interference terms with respect to the question, which and how many observables appear 
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FIG. 2: Diagrammatic representation of groups of longitudinal observables determining the interference terms of t-matrix 
elements for the helicity basis. The nomenclature for the groups and the corresponding observables are listed in Table IX. 



TABLE IX: Nomenclature for the diagrammatic representation of groups of longitudinal observables in Fig. 2 determining the 
interference terms of t-matrix elements for the helicity basis, where f[ M {X) is represented by X IM . 

Panel Line type Observables 

~) ^olid x{°, xf, yT, yT, xz 00 , xz 20 

(b) solid l 11 , l 21 , zf, zf, zf, zf, zz 11 , zz 21 

(c) solid l 22 , zf, zf, zz 22 

dashed xf, x\ 2 , y\ 2 , xz 22 

(d) solid xf, xf, xf, xf, yf, yf, xz 11 , xz 21 



in the representation of an interference term by observables. Because a closer inspection of the explicit expressions 
reveals, that in general the interference terms can be divided into disjunct subgroups which are determined by a 
subgroup of observables. In order to visualize this grouping we have devised in [6] a graphical representation. To this 
end one assembles the numbers "1" through "n" , where n denotes the number of i-matrix elements, by points on a 
circle and represents an interference term t*,tj by a straight line joining the points "j" and "j"'. Interference terms 
belonging to the same subgroup are then represented by the same type of lines. An example for the longitudinal 
matrix elements in the helicity basis is shown in Fig. 2. 

We will call a set of interference terms connected if they generate a pattern of connected lines so that any point 
belonging to one of the considered interference terms is connected to any other point of the set either directly or via k 
other intermediate points of that set. For example, in Fig. 2(a) one notes three disconnected lines, (b) and (d) contain 
two different groups of connected lines, whereas (c) contains one connected group of four lines. In such a connected 
set, any matrix element ty can be expressed in terms of any other matrix element tj of that set by one of the two 
forms 

*" . . . T iu-iiu-z !Wi t for k odd; 

V — \ T,.,, T,„,„ Ti, , Tj, „ Ta. „ , l lu 'J 



for k even, 



depending on whether the number k of intermediate points connecting j' with j via the points ji through jk is odd 
or even. The proof of these equations is simple and given in [6]. 

A special feature of the evolving geometric pattern is a closed loop, see e.g. Fig. 2(c). If such a closed loop has an 
even number of points then one finds from (107) for j' = j and k odd the following condition 

Tj k j k -iTj k _ 2 j k _ 3 ■ ■ ■ Tj 3 j 2 Tj 1 j = Tj 1 j 2 Tj 3 j 4 ■ ■ ■ Tj k 2 j k l Tj k j , (108) 

which means that in such a closed loop any interference term is completely determined by the other remaining inter- 
ference terms of that loop. This condition thus constitutes a complex relation between the participating observables 
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which allows one to eliminate two observables. On the other hand, a closed loop through an odd number of points 
(107) yields for / = j and k even 

Tjj k Tj k _ 1 j k _ 2 ' ' 'TjsjiTjij = TjjTj 1 j 2 Tj 3 j i ■ ■ ■Tj k3 j k2 Tj kl j k , (109) 

which again allows the elimination of two observables. It means furthermore that the modulus of each of the partici- 
pating t-matrix elements is completely determined, namely one has 

I J. |2 _ ± 3l3 ± J332 - L Jk-33k-4 ± ]k-l]k-2 rp 

1 1 j I _ rp rp T T 

Thus one may choose one matrix element of that loop as real and non-negative, fix its modulus according to (110) and 
then all other matrix elements of that loop are uniquely determined. As a side remark, we would like to point out, 
that the conditions in (108) and (109) constitute particular nonlinear relations between observables, which follow from 
the ones discussed in the Appendix B. In fact, these conditions are obtained by applying successively the condition 
in (B4) of this appendix to the left-hand sides of (108) and (109) yielding then the corresponding right-hand sides. 
As next step one has to choose from the total number of all interference terms t*,tj with j' > j which is ^n(n — 1) 
not counting t*,tj with f < j, because {t*,tj)* — t*jtj> - a set of n — 1 independent interference terms. Hereby we 
define a set of independent interference terms by the property that they generate a geometric pattern which does not 
contain any closed loop. From this definition follows that a set of n — 1 independent interference terms is represented 
by a pattern of n — 1 lines in such a fashion that (i) each of the n points is endpoint of at least one line, and (ii) each 
point is connected to all other points not necessarily in a direct manner but via intermediate points. It is obvious 
that in such a pattern no closed loops can be present, because one cannot construct from n — 1 lines a pattern which 
contains a closed loop and which still connects all n points. For such a set of n — 1 independent interference terms all 
matrix elements can be expressed by one arbitrarily chosen matrix element, say tj Q according to (107). In order to 
fix the remaining undetermined matrix element tj one has to choose one additional observable fo- From (105) one 
obtains in general an equation of the type 

fo = a + b\t JO \ 2 +c\t 30 \- 2 , (111) 

from which tj can be obtained, although not uniquely in general. The ideal situation would be such that one finds 
n — 1 independent interference terms each of them represented by only two observables. Because in this case one 
employs just 2n — 1 observables. On the other hand, analyzing the grouping of observables mentioned above, one 
will in general not find such a situation, either the number of observables for a set of n — 1 independent interference 
terms is larger than 2n — 2, or the grouping is such, that the choice of n — 1 independent interference terms involves 
observables which govern at least one additional interference term leading to one or several closed loops. However, in 
that case those loops lead to the elimination of superfluous observables. For example, considering Fig. 2, one notes 
that combination of the three lines of (a) with two connected lines of the connected groups in (b) or (d) results in 
a set of independent interference terms. According to Table IX, such a choice involves 14 structure functions, which 
means, that all six longitudinal matrix elements can be expressed in terms of 15 observables. In order to eliminate 
four of them, one can include the other group of two connected lines of (c) or (d), yielding two closed loops which 
allow the elimination. Thus, at the end one has a complete set of eleven observables. Further illustrative examples 
for such an analysis are discussed in some detail in [6]. 



(110) 
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IV. INGREDIENTS OF CALCULATION 

In this section we will review briefly the various ingredients which go into the calculation of the structure functions 
presented in this work. The structure functions are determined by hermitean quadratic forms in the matrix elements 
of the e.m. current operator between the initial and final states. As already mentioned in Sect. II, the calculation 
is done with respect to the cm. frame of the final hadronic NN state. Thus for a momentum transfer q the initial 
deuteron moves in this reference frame with a momentum —q, and one has to take into account the transformation 
to the lab frame as governed by the (3- factor of (37) . 

Structure functions are calculated within a nonrelativistic framework from the t- matrix defined in (27). Therefore, 
the wave functions are purely nonrelativistic and are obtained by solving the two-body Schrodinger equation with 
a realistic TV TV-potential for the bound as well as for the scattering states. Thus one of the principal ingredients of 
our calculation is a realistic potential model in order to generate the bound and scattering np-wave functions. 

The leading current contribution is provided by the nonrelativistic one-body nucleon current. As another important 
ingredient, we consider subnuclear degrees of freedom (d.o.f.) related to meson exchange currents (MEC) and isobar 
configurations (IC). Because of the increasing importance of relativistic effects with increasing energy and momentum 
transfer, we include also relativistic contributions of leading order beyond the nonrelativistic current. 

In view of the fact that some realistic NN-potentials are defined in r-space while others are in p-space, we 
have employed two separate corresponding codes. The r-space code is an outgrowth of the work of Fabian and 
Arenhovel [23] incorporating improvements and additions, particularly with respect to the most important rela- 
tivistic spin-orbit current and other leading order relativistic contributions. In this code, isobar configurations are 
treated in a perturbative approach [24] even though we had developed in the past a coupled channel r-space code [25] . 
However, in that code the A-propagator could not be treated as exactly as in a p-space code. 

The p-space code is described in detail in [26] . Besides inclusion of the N A-configuration as the most important 
isobar configuration in a coupled channel approach, it allows furthermore an exact treatment of all leading order 
relativistic contributions to one- and two-body currents as well as the Lorentz boost for a one-boson-exchange N N- 
interaction model. In fact, this code was developed for the Bonn p-space potentials (OBEPQ models) of [27] which 
are potentials of this type with short range cut-off form factors. 

A further difference between the two codes lies in their treatment of the electric multipoles. The r-space version 
incorporates the Siegert operators thereby insuring the inclusion of the dominant MEC implicitly in electric transi- 
tions. In the p-space version, mainly for historical reasons, the Siegert form of the electric multipoles is not used. In 
view of the fact that some potentials, e.g. the Bonn r-space, are derived from their corresponding p-space versions, 
one can estimate the inherent numerical differences between the two codes by using in this case the same potential 
model in both codes. 

The calculation of the i-matrix elements is based on an expansion of the final state into partial waves with total 
angular momentum j. The final state interaction (FSI) is taken into account by solving the corresponding scattering 
equation for a given partial wave \p, jmj). Then the electric and magnetic multipole transitions into this state are 
evaluated explicitly. For a given j the contributing multipoles are L = j — 1, j, j + 1. However, at some point one 
has to truncate the series at a maximal angular momentum j max for the explicit inclusion of FSI. On the other hand, 
the convergence of the partial wave expansion depends on the kinematics, in fact the convergence is quite slow for 
the quasi-free case. The solution to this dilemma is based on the fact that for higher partial waves the influence of 
FSI becomes increasingly unimportant with growing j so that for these the undistorted partial waves can be used 
instead. Thus we include all electric and magnetic multipoles up to a maximal multipolarity L max with consideration 
of FSI up to j max = L max + 1 and subtract the corresponding transitions without FSI and add finally the complete 
t- matrix obtained with a plane wave as final state, which we call the plane wave Born approximation (PWBA). This 
is described in detail in [23]. For all kinematics considered in this work L max = 4 was found to be sufficient, i.e. the 
final result did not change if we increased L max . As already mentioned, the higher partial waves are needed only in 
the vicinity of the quasi-free ridge because of the then slow convergence of the partial wave expansion. The formal 
expression for the t- matrix in PWBA is given in the Appendix F. We will now explain some details of the separate 
ingredients. 

A. AW-potential models 

As mentioned above, the two-body wave functions, needed for the calculation of the observables, are based on real- 
istic N A-potentials. In past work [1, 2, 4] we have employed a number of realistic potentials, such as Nijmegen [28], 
Paris [29], Bonn (r- and p-space versions) [27], and Argonne Vi 4 [30] potentials. In general, we found that the 
dependence of the observables on the choice of a realistic potential is rather moderate, in particular at low and medium 
energy and momentum transfers. For this reason we have chosen only one semi-modern potential, the Bonn p-space 
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model, for which a consistent meson exchange current has been constructed recently including all leading order rel- 
ativistic contributions [26] to the current operators, boost and internal dynamics as explained below in Sect. IV D. 
This potential, as most of the realistic iVA-potentials, is defined in purely nucleonic space without explicit A-d.o.f. 
Since, however, we will also present results, where such A-d.o.f. are treated explicitly in a coupled channel approach, 
we use in addition an interaction potential which has been constructed recently for this purpose [26] and which is 
based on the Bonn p-space potential. Furthermore, we consider as a prototype of recent high precision potentials the 
Argonne Vis potential [31] in the r-space code. 

B. One-body currents 

The one-nucleon current is derived from the nonrelativistic reduction of the Dirac current, retaining the leading 
order relativistic contributions. The internal nucleon structure is taken care of by including the free on-shell nucleon 
cm. form factors. In the p-space code we use the Dirac and Pauli form, Fi(Q 2 ) and F2(Q 2 ), respectively, where 
Q 2 = —q 2 . Explicit expressions are listed in the Appendix of [26]. The Sachs form with Ge(Q 2 ) for the electric and 
Gm(Q 2 ) for the magnetic form factor is used in the r-space code. It is obtained from the Dirac-Pauli form by the 
transformation 

Fi(Q 2 ) - t(Q 2 )(G e (Q 2 ) + ^G m (Q 2 )), (H2a) 

F 2 {Q 2 ) = t{Q 2 ){G e {Q 2 )-Gm{Q 2 )), (H2b) 

with t(Q 2 ) = (1 + Q 2 /4M 2 ) -1 and neglecting terms of higher order beyond the leading relativistic order. Since the 
difference between these two forms of the one-body current is of higher relativistic order, the different treatments in 
the two codes does not matter as long as the leading order relativistic contributions are included and as long as the 
kinematics stays within the limits of validity of the truncated expansion keeping only the terms of leading relativistic 
order [32]. 

A variety of form factor parametrizations is available [33-37]. The largest uncertainty exists for Geu- Indeed, one 
of the early motivations for investigating deuteron electrodisintegration with polarization degrees of freedom was that 
it could provide a nearly model independent method of determining the neutron electric form factor Geu [38-40]. 
Since in these studies the sensitivity of various observables with respect to Geu has been investigated extensively, 
we have employed here only one form factor model, namely the dipole model including a nonvanishing electric form 
factor of the neutron in the Galster parametrization [33] (with p=5.6). 

In principle, one should consider also off-shell effects in the one-body current in view of the fact, that the nucleons 
are not free but subject to the hadronic interaction. In such a situation, the form factors would acquire an additional 
dependence on the initial and final squared four-momenta of the nucleons and, moreover, additional currents with 
more off-shell form factors would appear. That such effects potentially may be non-negligible has been shown recently 
for deuteron photodisintcgration [41]. In that study off-shell effects were evaluated using a simple, but not very 
realistic pion cloud model for the nucleon structure. At present, however, no realistic treatment of such off-shell form 
factors exists, and thus we neglect them here. 

C. Meson exchange currents 

An important property of realistic AA-interaction models is that they induce two-body meson exchange currents 
(MEC). Indeed, any isospin and/or momentum dependence of an TV TV-potential requires on a formal basis the existence 
of an interaction current in order to satisfy the continuity equation. The physics underlying such MECs is related to 
the coupling of the hadronic interaction diagrams to the e.m. field. In the case that the AW-potential is explicitly 
derived from a meson-exchange model, this connection is obvious and one obtains straightforwardly the associated 
nonrelativistic MEC, of which the 7T-MEC is the most important one, as well as the leading order relativistic two-body 
charge and current contributions. 

However, for potential models, which in their medium range part use a phenomcnological parametrization, the 
construction of a proper exchange current for the isospin dependent potential part is not unique, because the connection 
to the underlying physical process is obscured. For such cases, a recipe has been developed in the past independently 
by Riska [42] and by Buchmann et al. [43], which is inspired by the genuine meson exchange models. This recipe 
is based on the observation that the spin-isospin dependent central and tensor parts of a given N iV-potential can be 
split into a pion- and rho-exchange-like potential for which the corresponding meson exchange currents are known. 
While the approach of RiSKA is based on the momentum space representation of the potential and thus can be applied 
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to any phenomenological potential, the method of Buchmann et al. was conceived for an r-space representation of 
the potential as a superposition of appropriate Yukawa functions, and thus its application appeared to be limited to 
such type of potentials. But recently, this approach was extended to potentials with a more general radial behaviour 
by applying a Laplace transform [44]. 

This phenomenological method works for the nonrelativistic MEC reasonably well, but one should be aware, that its 
construction contains some inherent arbitrariness, because the it- and p-MEC contain purely transverse pieces, which 
are not constrained by current conservation and thus can be modified arbitrarily without destroying the consistency. 
An example for such a modification is given in [43]. Thus, extending this recipe to the construction of relativistic 
MEC contributions in a corresponding manner has to be considered as purely heuristic. Fortunately, an important 
part of MEC can be incorporated model independently by the use of Siegert operators for the transverse electric 
multipoles [45, 46]. In fact, in this way the major MEC contribution to electric transitions is consistently incorporated 
implicitly, at least for low and medium energy and momentum transfers. In fact, the results in our work [1, 2, 4, 5] 
labeled as the normal part (N) contain MEC contributions via the electric Siegert operators implicitly, but otherwise 
no explicit MEC in either the electric or magnetic transitions. However, as already mentioned, in the p-space code, 
based on the work of Ritz ct al. [26], no Siegert operators are used. In this case the results labeled "N" do not 
contain any MEC implicitly. With respect to the isobar configurations, which will be discussed in the next section, 
we also include those isobar-MEC which are induced by the transition potentials. 

A last remark concerns the question of e.m. form factors for the MEC. Here we take the heuristic approach 
multiplying the isoscalar and isovector pieces by the appropriate isoscalar and isovector nucleon form factors. In this 
way, current conservation holds also in the presence of such form factors. 

D. Isobar contributions 

Isobar d.o.f., describing phenomcnologically internal nucleon d.o.f., can be incorporated either in the form of effective 
nonlocal two-body operators, describing intermediate excitations of one or two isobars, or by allowing explicit isobar 
configurations (IC) in the nuclear wave functions, where one or several nucleons are replaced by an isobar, and with 
appropriate strong and cm. operators [47]. In the present work the latter approach is used by admitting isobars as 
explicit constituents in the two-body system. This allows one also to handle in a natural manner the real excitation 
of a A(1232)-resonance for high enough energy transfers above pion production threshold. It furthermore avoids the 
often applied static approximation of the effective MEC induced by the intermediate excitation of an isobar, which 
has very limited value only. In the present work, we have included the AA, AA(1440), and AA configurations in the 
r-space code while in the p-space code only the AA configuration is included as the most important one for energy 
transfers up to about 400 MeV. 

The corresponding wave function components, called isobar configurations, are obtained either in a perturbative 
treatment or in a coupled channel approach. In the perturbative approach, an isobar configuration, consisting, for 
example of a AA configuration, is generated by just one AA-collision via a transition potential AA — ► AA for 
which a simple one-boson exchange model is used [47] . On the other hand, in a coupled channel approach one has to 
renormalize the original AA-interaction as mentioned above, because being fitted in pure AA-space to experimental 
scattering data, it contains implicitly already the effect of such intermediate configurations, e.g. AA. This means 
that in principle one would need to redo the fit of the potential parameters if such isobar configurations are included 
explicitly In order to avoid such involved work, a reliable box-subtraction method, first proposed by Green and 
Sainio [48], is available and was applied in [25, 26]. In most cases, however, we use the simpler perturbative 
calculation [47]. Only in case of the AA-configuration do we also consider a coupled channel calculation for the Bonn 
p-space potential [26] as mentioned above, because in momentum space the A propagation can be treated in a more 
exact manner compared to an r-space calculation [25] . Finally, one- and two-body current operators involving isobars 
have to be considered. In view of the fact, that these currents are less well known, we neglect relativistic terms and 
restrict ourselves to their nonrelativistic expressions which are given in [47] for the one-body terms and in [49] for the 
MEC contributions, again with appropriate e.m. form factors. 

E. Relativistic contributions 

Relativistic contributions arise from three sources. These are (i) the internal relativistic dynamics in the rest 
frame of the nucleus, (ii) the boost of the intrinsic wave function from the nuclear rest frame, here the lab frame, 
to a moving one because of the nonvanishing momentum transfer, and finally (iii) relativistic contributions to the 
interaction operators, here the current operators. In the present approach we resort to a p/M-expansion, retaining 
only the leading order relativistic terms beyond the nonrelativistic limit. The boost of a wave function is described 
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by a unitary operator for which the p/M-expansion yields in leading order two separate contributions, a purely 
kinematic part, which can be interpreted as the effect of Lorentz contraction and Thomas- Wigner spin rotation, 
and a potential dependent part, which only is present for pseudoscalar meson exchange [50]. Relativistic contributions 
to the one-body current have been discussed already in Sect. IV B. With respect to MEC, consistent treatments are 
available in [51-54] which are based on a meson-theoretical one-boson-exchange potential as nuclear interaction. The 
construction is more questionable for semi-phenomcnological potentials like the Argonne Vis, as pointed out in 
Sect. IV C. Even a consistent nonrelativistic MEC is ambiguous for such potentials [42]. Although one can proceed 
also for the relativistic MEC in analogy to purely one-boson-exchange models, one should be aware of the inherent 
ambiguities of such an approach. 

A consistent treatment of all three types of contributions is given in the work of Goller and Arenhovel [52] for 
a pure one-pion-exchange model and by Tamura et al. [53] for a more general potential type. The work of [52] has 
been generalized in [26] to a consistent leading order relativistic treatment for the p-space Bonn potentials, and our 
p-space code is based on this work. In view of the problems associated with the construction of a consistent MEC 
for a semi-phenomcnological potential, we include in the r-spacc code as relativistic current contribution only the 
one-body part, containing the most important spin-orbit current, and the kinematic boost as described in [55]. 



In the following discussion we will use as shorthand "FSI" for final state interaction, "PWBA" for nonrelativistic 
one-body current with plane wave final state, "RPWBA" for relativistic one-body current with plane wave final state, 
"N" for the nonrelativistic normal theory, i.e. without explicit meson exchange currents and isobar configurations. 
This means in the case that Siegert operators are used, that a part of MEC is implicitly included in the electric 
multipoles for "N" as is the case for the r-space code, whereas for the g-space code no MEC contributions in "N" 
appear. Furthermore, "MEC" and "IC" stand for the contributions of explicit meson exchange currents and isobar 
configurations, respectively, "RC" for the inclusion of relativistic contributions and "T=N+MEC+IC+RC" for the 
complete calculation. 



The inclusive reaction d(e, e')np for unpolarized beam and target is governed by a longitudinal form factor Fl and a 
transverse one Ft, whereas eight additional form factors appear if one allows for beam and target polarization. These 
form factors depend on two variables, for which we choose the final state cm. excitation energy E np (see Eq. (1)) and 
the squared cm. three-momentum transfer q 2 (from now on we use the notation q = <f c m ), i.e. F L / T = F L / T (E np , q 2 ). 
All results presented in Figs. 3 through 11 arc obtained using the Argonne Vis potential. In order to give an overview, 
the upper panels of Fig. 3 show Fl and Ft, calculated with the complete theory, in combined surface and contour 
plots over the E np - q 2 -planc for E np = — 300 MeV and q 2 = — 25 fm~ 2 . For a more detailed view of the form 
factors in the near threshold region the lower panels display them in a smaller part of the E np - q 2 -plane, namely for 
E np = - 20 MeV and q 2 = - 2 fm" 2 . 

One readily notes that the quasi-free ridge along E np /MeV~ 10<7 2 /fm~ 2 is the dominant feature of these two form 
factors, where the quasi-free kinematic is defined by the requirement that the virtual exchanged photon is absorbed 
by only one nucleon, which is emitted in the forward direction, with energy and momentum transfer such that the 
spectator nucleon remains at rest in the lab system. This yields the condition 
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which gives for the invariant mass according to (4a) using (9b) 
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or for the final state cm. excitation energy 
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FIG. 3: Surface and contour plots of longitudinal and transverse form factors as function of E np and q 2 for E np = — 300 MeV 
and q 2 — - 25 fm~ 2 using the Argonne Via potential. 




E n [MeV] E n [MeV] 



FIG. 4: Upper panels: Longitudinal and transverse form factors along the quasi-free ridge as function of E np calcu- 
lated using the Argonne Vis potential. Notation: dotted: PWBA; long dashed: RPWBA; dash-dot: normal (N), i.e. 
nonrelativistic approach with FSI included; short dashed: nonrelativistic MEC and IC included (N+MEC+IC); solid: 
complete calculation (T=N+MEC+IC+RC). Lower panels: Form factor ratios R L/T = F i/T (N)/F L/T (PWBA) (dot- 
ted), R L/T = Fl /t (T) / Fl it (RPWBA) (long dashed) R L/T = F L/T (N+MEC+IC) /F L/T (N) (dash-dot), and R L/T = 
F L/T (T)/F i/T (N+MEC+IC) (short dashed). 

The latter relation gives the already mentioned rule of thumb E np /MeV~ 10 g 2 /fm -2 , a straight line in the E np - q 2 - 
plane as is obvious in Fig. 3. 

The behaviour of the form factors along the quasi-free ridge is displayed in the upper two panels of Fig. 4 while 
the lower two panels show form factor ratios with respect to the various ingredients and interaction effects. The 
longitudinal form factor rises steeply, reaches its maximum at quite low E np around 5 MeV, and falls off very 
rapidly with increasing E np . On the other hand, Ft rises considerably slower and reaches its maximum only around 
E n p — 30 MeV. Also the fall off is much slower compared to Fl ■ One furthermore notes that for Fl the effect of MEC 
and IC is unimportant while for Ft these are still sizeable above E np = 20 MeV and of the order of several percent 
decreasing slowly with growing E np as shown by the dash-dotted curves in the lower panels. Also the influence of 
FSI becomes quite unimportant above E np rts 30 MeV as seen in the upper panels by comparing the dotted curves 
(nonrelativistic PWBA) with the short dashed curves (N+MEC+IC). Quantitatively one finds from the ratios in the 
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FIG. 5: Surface and contour plots over the E np -q 2 -plane of ratios of longitudinal form factors for different contributions: (a) 
influence of FSI, Fl (PWB A) / Fl (N) , for low E np = - 20 MeV and in (b) for E np = 20 - 300 MeV; (c) influence of IC, 
F L (MECIC)/F L (MEC), for low q 2 = - 2 far 2 and in (d) for q 2 = 2 - 25 fai" 2 ; (e) influence of RC, F L (T)/F L (MECIC). 



lower panels (dotted curves) that above E np w 100 MeV an almost constant difference of a few percent remains. The 
analogous ratios with relativistic contributions included, i.e. F L / T (T) /F L i T (RPWB A) show a similar behaviour (long 
dashed curves). The reason for the large overestimation of Fl in PWB A at low E np has its origin in the fact that 
the final state plane wave is not orthogonal to the deuteron bound state so that the charge monopole transition is 
not suppressed near threshold. On the other hand, Ft is strongly underestimated in PWBA because of the absence 
of the resonance in the So state. The only notable effect arises from relativistic contributions leading for both form 
factors to a sizeable reduction which increases almost linearly with E np . However, these RC are quite well accounted 
for in the relativistic RPWBA as is demonstrated by the little difference between the complete calculation and the 
relativistic RPWBA their ratios approaching one with increasing E np (long dashed curves). 

We now will turn to the near threshold behaviour shown in the lower panels of Fig. 3. One readily notes that for 
E np — * along g 2 =const. Fl runs first through a broad maximum and then decreases rapidly to zero while Ft rises 
dramatically resulting in a very sharp peak right above threshold. The peak height grows first with increasing q 2 , 
reaches its maximum around q 2 = 0.5 fm~ 2 and then falls off. The rapid decrease of Fl and the sharp peak of Ft 
near E np w is a consequence of the fact, that close to break-up threshold the 1 So-scattering state, the so-called 
anti-bound state, dominates the final state into which the Coulomb monopole transition is forbidden while one has 
a very strong magnetic dipole isovector transition which is further enhanced by MEC and IC contributions. Above 
E np 4 MeV and q 2 w 0.4 fm~ 2 both the surface plot as well as the contour lines exhibit clearly the onset of the 
quasi-free ridge in both form factors. 

The relative influence of the various interaction effects are shown as ratios in Fig. 5 for the longitudinal form factor. 
Because of the strong influence of FSI on Fl near threshold, we show for this form factor the ratio Fl (PWBA)/_FY, (N) 
separately for the region E np < 20 MeV in Fig. 5 (a), exhibiting the already noted very strong effect of FSI, resulting 
in a strong increase of this ratio for q 2 — ► along E np =const., whereas Fl(P~WBA) / Fl(N)—* for E np — > along 
q 2 =const. For the remaining part (20 < E np /MeV < 300) this ratio is shown in Fig. 5 (b). The influence of FSI is 
minimal on top of the quasi-free ridge while it leads to an increase above this ridge, i.e. for E np /MeV > 10q 2 /im~ 2 , 
and to a decrease below. MEC have almost no effect on Fl and thus are not shown here, because the dominant 
nonrelativistic 7T-MEC does not contribute to the charge density. The effect of isobar configurations on Fl is shown 
in Fig. 5 (c) for q 2 < 2 fm -2 and in (d) for q 2 > 2 fm -2 . Above the quasi-free ridge Fl is reduced by IC. In particular 
for low q 2 close to zero and E np approaching the A-excitation region, Fl decreases drastically by IC to about 10 
percent. Going below the quasi-free ridge, one again notes a reduction but of smaller size which diminishes when 
approaching small E np . Furthermore, RC, shown in Fig. 5 (e), exhibit an interesting behaviour: along the quasi-free 
ridge one finds a distinctive valley describing the increasing reduction by RC with increasing E np or q 2 as was already 
apparent in Fig. 4 (upper right panel). Away from the quasi- free ridge the influence of RC diminishes first on both 
sides to almost zero and increases then again when approaching the regions of higher E np and lower q 2 or vice versa. 
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FIG. 6: Surface and contour plots over the F„ p -g 2 -plane of ratios of transverse form factors for different contributions: (a) 
influence of FSI, Ft (P WB A) /Ft (N) ; (b) influence of MEC, Ft (MEC) /Ft (N) for low E np = - 20 MeV and in (c) for E np = 
20-300 MeV; (d) influence of IC, F T (MECIC)/F T (MEC); (e) influence of RC, F T (T)/F T (MECIC) for low E np =0-20 MeV 
and in (f) for E np = 20 - 300 MeV. 

The corresponding ratios for Ft are displayed in Fig. 6. For this form factor the effect of FSI is much more 
pronounced off the quasi- free ridge than in Ft as is shown by the ratio Ft (P WB A) /Ft (N) in part (a). One finds a 
very strong decrease on both sides, i.e. at low E np with increasing q 2 as well as at low q 2 but increasing E np , which 
means a strong enhancement by FSI. For the display of MEC effects we show in part (b) the region of low E np and in 
part (c) the remaining region. As already mentioned, along the quasi-free ridge one finds little influence. But going 
away from this ridge, MEC lead to a sizeable increase as shown in part (c), especially strong close to threshold (see 
part (b)) up to about q 2 = 15 fm -2 . For higher momentum transfers MEC result in a reduction. This behaviour 
is well known and in agreement with experimental data [56, 57]. Isobar effects displayed in Fig. 6 (d) become quite 
pronounced only in the region of A excitation near E np = 260 MeV and for not too high momentum transfers. Finally, 
relativistic contributions lead for E np near threshold and q 2 around 15 — 20 fm -2 to quite a significant increase as 
shown in Fig. 6 (e). For the remaining region in Fig. 6 (f) one notes again a significant reduction along the quasi- free 
ridge as in Ft in the form of a pronounced valley. 

A survey for the additional form factors for polarized beam and target, calculated with the complete theory, is 
shown in Fig. 7 for F 20 , F 20 , F T °, and F^ T 2 and in Fig. 8 for the iF-interference form factors as surface and 
contour plots over a smaller portion of the E np - g 2 -plane, i.e. E np = — 160 MeV and q 2 = — 16 fm~ 2 . The largest 
polarization form factors are F T °, F^ T 2 , F r °, and F'It 1 which are of the same order of magnitude as the unpolarized 
form factors. An order of magnitude smaller are F|° and F 2 ^ 1 . The remaining two, F^ 1 and F'^f 1 , are three orders 
of magnitude smaller although they increase slightly in size along the quasi-free ridge. 

The influence of FSI, MEC, IC and RC cannot be displayed as ratios because these form factors have zeros. For 
this reason we show only for the larger form factors the most important influences by plotting differences for FSI 
(N-PWBA) in Fig. 9, for combined contribution from MEC and IC (MECIC-N) in Fig. 10, and for RC (T-MECIC) in 
Fig. 11. Final state interaction is very important for all polarization form factors shown in Fig. 9 in the region of low 
E n p = — 40 MeV and low q 2 = — 4 fm~ 2 . MEC and IC effects displayed in Fig. 10 arc substantial in the transverse 
form factors F^P, F TT , and F T °, mostly below the quasi- free ridge, whereas they are of minor importance in the 
interference form factors F 2 ^ 1 and F'^f 1 , the largest effect being near the quasi-free ridge. Relativistic contributions 
show up in Fig. 11 along the quasi- free ridge in the primed form factors F T ° and F' LT ~ as well as in F L ^, . 

Now we will turn to a comparison with experimental data. Most extensively discussed has been the inclusive 
electrodisintegration near threshold, averaged over E np = through 3-5 MeV, at higher momentum transfers (see [56, 
57]) and we will not repeat the discussion here. Instead we will compare the theory to the experimental data for 
Ft and Ft as obtained by a Rosenbluth separation in the near threshold region for various momentum transfers by 
Simon et al. [58]. The only comparison to theory reported in that work was for a nonrelativistic treatment employing 
the Paris potential and 7T-MEC, and a satisfactory agreement was found. In the light of the considerable progress 
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FIG. 7: Surface and contour plots of polarization form factors F^ , F^°, F^°, and F^, T 2 for polarized beam and target as 
function of E np and q 2 for E np = - 160 MeV and q 2 = - 16 fin" 2 . 




FIG. 8: Surface and contour plots of polarization LT-interference form factors for polarized beam and target as function of 
E np and q 2 for E np = 0-160 MeV and q 2 = - 16 fm -2 . 

which theory has achieved since then, it appears timely to compare the modern approaches with those data. Thus 
we confront in Fig. 12 three data sets of [58] corresponding to the squared momentum transfers q 2 w 0.6, 1.5, and 
3.8 fm -2 for E np = — 9 MeV with various interaction effects, i.e. normal nonrelativistic approach (N), with inclusion 
of meson exchange currents and isobar configurations (N+MEC+IC), and the complete theory including relativistic 
contributions (T=N+MEC+IC+RC) for the BONN-Qb potential. 

As expected, for Fl the interaction effects from MEC and IC are very small though not completely negligible. 
Relativistic effects show some influence resulting in a small reduction above E np w 5 MeV which increases slightly in 
size with E np . The agreement of the full theory with the data for Fl is quite satisfactory. One should keep in mind, 
that no open parameter has been fit in this comparison. 

The transverse form factor Ft shows right above threshold the already noted distinctive peak arising from the 
Ml-transition into the anti-bound 1 So-state. Here the influence of MEC and IC is quite significant, especially in the 
peak region, and it grows sizeably with increasing squared momentum transfer. On the other hand, RC lead only to 
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FIG. 9: Influence of FSI: Surface and contour plots of differences of polarization form factors calculated in normal nonrelativistic 
theory (N) and in PWBA as function of E np and q 2 for E np = 0-140 MeV and q 2 = - 14 fm~ 2 . 




FIG. 10: Influence of MEC and IC: Surface and contour plots of differences of polarization form factors calculated with MEC 
and IC and in normal nonrelativistic theory as function of E np and q 2 for E np = — 140 MeV and q 2 = — 14 fm~ 2 . 

a slight reduction for q 2 = 3.8 fm~ 2 while for the other two cases they show no effect at all. The agreement with the 
data is in general quite good except for the peak where the theory lies above the data. It is very likely, however, that 
the experimental resolution was insufficient to resolve this very sharp peak. The potential model dependence is shown 
in Fig. 13 using the BONN-Qb, the two Argonne potentials V14 and Vis, and the Paris potential. The variation 
of the predicted form factors by these potentials is much smaller in size than the size of interaction effects. The only 
exception is the prediction of the peak height for the Paris potential, which gives a lower value. The reason for this 
lies in the Paris potential's prediction of too small a value for the np-scattering length. Finally, in order to give a 
more detailed and quantitative comparison with experiment we show in Fig. 14 the ratios of the data and the various 
theoretical model predictions to the results obtained with the Bonn-QB potential. It is obvious that the variation 
with the potential model is substantially smaller than the experimental errors. Certainly, much more precise data, in 
particular at the threshold peak, are needed in order to put the theory to a more critical test. 
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FIG. 11: Influence of RC: Surface and contour plots of differences of polarization form factors calculated with MEC and IC 
and for complete theory as function of E np and q 2 for E np — — 140 MeV and q 2 — — 14 fm~ 2 . 





FIG. 12: Longitudinal and transverse form factors with var- 
ious interaction effects. Notation: dash-dot: N; dashed: 
N+MEC+IC; solid: N+MEC+IC+RC. Experimental data 
from Simon et al. [58]. 



FIG. 13: Longitudinal and transverse form factors for vari- 
ous potential models. Notation: dashed: Paris; long dash- 
dot: Argonne V14; short dash-dot: Argonne Vis; solid: 
BONN-Qb. Experimental data from from Simon et al. [58]. 



Another set of data for the inclusive reaction at higher excitation energy E np and higher momentum transfers is 
provided by Quinn et al. [59]. We show first a comparison between theory and experiment for some inclusive cross 
sections corresponding to several different kinematics reported by Quinn et al. which exhibit nicely the quasi-free 
peak, namely at E np w 20 MeV in Fig. 15 and at E np w 70 MeV in Fig. 16. The upper panels refer to forward angles 
(60°) and the lower ones to backward angles (134.5°). In the left panels the predictions of the nonrelativistic normal 
theory obtained with the BONN-Qb and Argonne Vis potentials is compared to the data and in the right panels 
the complete theory. The normal theory results in a slight potential dependence as is apparent in the quasi-free peak 
with the results for the BONN-Qb potential slightly higher and more pronounced for the 60°-data. However, this 
potential model dependence is very much reduced for the complete theory. The most significant improvement of the 
full calculation is seen in the near threshold region at backward scattering angles (see lower panels). In general the 
agreement is quite good. The underestimation above u = 200 MeV in the upper panels of Fig. 16 has its origin in the 
absence of pion production contributions in the theory which become significant in this region. 

The fact, that the model dependence of the nonrelativistic normal theory is stronger for forward scattering angles, 
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FIG. 14: Ratios of the longitudinal and transverse form factors for various potential models with respect to the BONN-Qb 
potential. Notation as in Fig. 13. Experimental data from from Simon et al. [58]. 
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FIG. 15: Inclusive cross sections d(e,e). Upper panels 
for £} ab = 292.8 MeV and 6 l f> = 60°. Lower panels for 

E Ub = 174 3 Me y and e \ah = 134 go Exp from Q UINN et 

al. [59]. Notation of curves: Left panels for nonrelativistic 
normal theory (N) and right panels for complete theory (T) 
for BONN-Qb (dashed) and Argonne Vis (solid) potentials. 



FIG. 16: Inclusive cross sections d(e, e). Upper panels 
for E\ ah = 596.8 MeV and 6^ ab = 60°. Lower panels for 
M ab = 367.7 MeV and 6 l f° = 134.5°. Exp. from Quinn et 
al. [59]. Notation of curves: Left panels for nonrelativistic 
normal theory (N) and right panels for complete theory (T) 
for BONN-Qb (dashed) and Argonne Vis (solid) potentials. 



points to a stronger model dependence of Fl compared to Ft- Indeed, this is confirmed by the comparison of the 
theory with the experimentally determined longitudinal and transverse form factors in Figs. 17 and 18. However, this 
model dependence disappears almost completely for the full theory. 

One should keep in mind that a direct comparison with the response functions R L / T reported in [59] is not possible 
because of different definitions. However, the relations of Rl/t to our form factors is easily obtained by comparing 
the formal expressions for the inclusive cross section. In Ref. [59] the response functions Rl/t & re defined by 

^F = <7Mott (^ + (r? + 2 )iiT )' (116) 
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FIG. 17: 



Longitudinal form factors. Upper panels for 
278.5 MeV and 9^ h = 134.5° and lower panels for 
134.5°. Exp. from Quinn et 
al. [59]. Notation of curves: Left panels for nonrelativistic 
normal theory (N) and right panels for complete theory (T) 
for BONN-Qb (dashed) and Argonne Vis (solid) potentials. 



Ef h = 174.3 MeV and 6»i ab 



FIG. 18: Transverse form factors. Upper panels for E l ± h = 
278.5 MeV and 6 l f° = 134.5°, middle panels for E l f = 
174.3 MeV and O 1 ^ = 134.5°, and lower panels for E[* h = 
367.7 MeV and 9^ = 134.5°. Exp. from Quinn et al. [59]. 
Notation of curves: Left panels for nonrelativistic normal 
theory (N) and right panels for complete theory (T) for 
BONN-Qb (dashed) and Argonne Vis (solid) potentials. 



whereas we use according to (78) in conjunction with (26) 

6c(fc> ab , k l * h )(p L F L + PT F T ) 

This yields finally 

a 2 

p L = 2 7t 2 a jjt— Rl an d F T = 2it 2 aR T . (118) 

For Fl the normal nonrelativistic theory reveals a sizeable potential model dependence according to the left panels of 
Fig. 17 which, however, is strongly reduced for the complete theory leading to a satisfactory agreement with the data. 
On the other hand Ft in Fig. 18 exhibits much less sensitivity to the potential model for the normal nonrelativistic 
theory The improvement by the full calculation in the near threshold region with respect to the data is clearly seen 
in the middle panels. In the lower panel the absence of pion production in the theory is again responsible for the 
underestimation above w = 200 MeV. 
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B. Exclusive Observables 



Compared to the eight inclusive observables below pion threshold one has 324 independent structure functions 
for the exclusive case. This greater variety coupled with the angular dependence of each observable allows a more 
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FIG. 19: Survey on the four unpolarized structure functions along the quasi-free ridge, calculated for the BONN-Qb potential. 
The top left inset "[(-n) fm]" indicates the unit [10~ n fm] for the structure function and the top right inset ll [E np /q 2 \\ where 
E np in [MeV] and q 2 in [fm -2 ], indicates the kinematic sector. Notation of the curves: N (short dashed), N+MEC (dash-dotted), 
N+MEC+IC (long dashed), total=N+MEC+IC+RC (solid). 



detailed analysis of the reaction under study. For example, the differential cross section alone provides 41 structure 
functions if one allows for beam and target polarizations (see Table III). Polarization analysis of one or both outgoing 
nucleons yields an even larger number of structure functions. However, as has been discussed in detail in Sect. IIIB 2, 
in principle a set of 35 independent observables should suffice for a complete determination of all reaction matrix 
elements. In practice, however, some of these observables will be very difficult to measure experimentally with the 
required accuracy. 

For this reason we have decided to discuss here only those structure functions which are relatively easily accessible, 
i.e. which require the measurement of only one asymmetry (see Table X of Appendix D). Our discussion will focus 
on the sensitivity of these structure functions to various interaction effects. 

To begin we will briefly give a survey on the four structure functions of the unpolarized differential cross sec- 
tion. Since the structure functions depend on E np and q 2 in addition to Q c m -^ we have chosen to represent the 
various kinematical regions by a grid in the E np — g 2 -plane, denned by E np = 10, 50, 100, 175, 250 MeV and 
q 2 = 1, 5, 10, 17.5, 25 fm~ 2 . Only for E np = 250 MeV the lowest g 2 -value was taken as 1.5 fm~ 2 because of 
the photon line according to (7). Fig. 19 shows the four structure functions along the quasi-free ridge. The lon- 
gitudinal structure function fi, shows a pronounced peak in the forward direction essentially caused by the charge 
interaction of the virtual photon with a proton which is emitted preferentially along q. There is no corresponding 
peak at 180° for the neutron because of its very small electric form factor. For this reason we have restricted the 
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FIG. 20: Survey on longitudinal structure function _/x off the quasi-free ridge, calculated for the BONN-Qb potential. Notation 
of the curves as in Fig. 19. 



angular range to the forward direction. With increasing momentum transfer, the width of this forward peak decreases 
markedly (one should note that the angular range differs for the various cases). Similarly, along the quasi-free ridge 
Jt exhibits forward and backward peaks with decreasing width for growing q. The forward peak arises from the e.m. 
interaction with the dominant spin current of the proton at 9 = 0° while the backward peak at 9 = 180° (again not 
shown) is similar in structure and corresponds to neutron emission along q. The ratio of the forward to backward 
peaks is essentially given by the square of the ratio of the proton to neutron magnetic moments i.e. (/j, p /^ n ) 2 ~ 2. 
As was already pointed out in the discussion of the form factors, along the quasi-free ridge most of the interaction 
effects arc marginal, except for relativistic contributions which can result in a decrease of the forward and, in the case 
of /t, backward peaks growing sizeably with increasing momentum transfer. The decrease is particularly significant 
in Jl amounting, for example, at q 2 = 10 fm~ 2 to about 40 % and at q 2 = 25 fm -2 to even 70 %. In fx the effect is 
smaller, roughly by a factor two. 

The interference structure function Jlt, which has to vanish for 9 = 0° and 180°, exhibits a negative forward 
peak of broader distribution than the diagonal structure functions and is almost an order of magnitude lower. A 
much smaller peak appears for neutron emission close to 180° with opposite sign which is not shown in Fig. 19. The 
influence of MEC and IC is very tiny and even RC show up only slightly at intermediate E np — 50 to 100 MeV. 
The interference structure function Jtt is even smaller, about two orders of magnitude compared to Jlt- Interaction 
effects from MEC and IC play only a slightly more significant role. However, RC result in a distinct and sizeable peak 
at backward angles. 

In Figs. 20 through 23 we show the four structure functions for kinematical settings off the quasi-free ridge, i.e. six 
settings for E np /MeV> 10g 2 /fm -2 in the upper left panels and six settings for _E„ p /MeV< 10<7 2 /fm~ 2 in the lower 
right panels. Away from the quasi-free ridge, the magnitudes of Jl and drop rapidly. For shown in Fig. 20, 
the width of the forward peak increases. On the lower right side only RC show quite sizeable influences because of 
the absence of nonrelativistic MEC contributions and the neglect of a very small charge excitation of the A - for this 
reason the short and long dashed and the dash-dot curves coincide in Fig. 20 -, whereas on the upper left side, i.e. 
for _E„p/MeV> 10g 2 /fm~ 2 , also IC become increasingly important when approaching the A-region while keeping the 
momentum transfer small. The influence of IC arises through the change of the normal wave function component via 
the dynamic coupling of the NN- and N A-channels. 

For the transverse structure function in Fig. 21, which is slightly larger in size than the various interaction 
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FIG. 21: Survey on transverse structure function fr off the quasi-free ridge, calculated for the BONN-Qb potential. Notation 
of the curves as in Fig. 19. 
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FIG. 22: Survey on interference structure function /lt off the quasi-free ridge, calculated for the BONN-Qb potential. Notation 
of the curves as in Fig. 19. 
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FIG. 24: Influence of IC on various structure functions calculated with the BONN-Qb potential. Notation: N (short dashed), 
N+MEC (dash-dotted), N+MEC+IC (long dashed), T=N+MEC+IC+RC (solid). 



effects become much more pronounced off the quasi-free ridge than for _/x because of the presence of nonrelativistic 
MEC and strong transverse excitation of the A. Indeed, at low energy MEC provide the largest interaction effect 
followed by RC while IC remain small. That changes for energies in the region of A-excitation above the quasi- 
free ridge, where naturally the IC contributions dominate, although MEC are also sizeable and even RC cannot be 
neglected, in particular in the forward and backward regions. 

The interference structure function Jlt in Fig. 22 is comparable in size to /l, however, considerably more sensitive 
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FIG. 25: Upper panels: Influence of MEC on various structure functions calculated with the BONN-Qb potential. Notation as 
in Fig. 24. Lower panels: Potential model dependence of the total result for the same structure functions. Notation: Bonn-R 
(short dashed), Paris (dash-dotted), Argonne Vi 8 (long dashed), BONN-Qb (solid). 
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FIG. 26: Upper panels: Influence of RC on various structure functions calculated with the BONN-Qb potential. Lower panels: 
Potential model dependence of the total result for the same structure functions. Notation as in Fig. 25. 



to interaction effects. It shows a negative peak in forward direction, a maximum around 90° and another but much 
less pronounced negative peak at backward angles. RC produce by far the strongest effect, deepening the forward 
negative peak, enhancing the maximum and weakening or even washing out completely the negative backward peak 
so that a small positive peak results for some kinematics. 

The other interference structure function Jtt, shown in Fig. 23 is of comparable size to /lt off the quasi- free ridge. 
This is in contrast to their behavior along the quasi-free ridge where Jtt is about two orders of magnitude smaller than 
Jlt as depicted in Fig. 19. Below the quasi-free ridge it shows a distinct minimum around 110° which is decreased 
by MEC but then deepened slightly by IC and, but more strongly, by RC. Above the quasi-free ridge the influences 
of interaction effects show in a certain sense an opposite behaviour. Here MEC deepen the minimum, lying more 
at forward angles, but produce a sign change in the backward direction. This effect is drastically counterbalanced 
at forward angles by IC and amplified at larger angles. The additional RC lead then finally to a smaller reduction, 
deepening the small forward minimum, which exists for higher E np and q 2 . 
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FIG. 27: Influence of IC on asymmetries for kinematic settings for the separation of the various structure functions in Fig. 24 
calculated with the BONN-Qb potential. Notation: N (short dashed), N+MEC (dash-dotted), N+MEC+IC (long dashed), 
T=N+MEC+IC+RC (solid). The upper left inset lists the various angles according to [8 e , <j>, 9d, 4>d]- 



Now we turn to the question of which structure functions exhibit the most significant sensitivity with respect to 
the various interaction effects. As already mentioned we take into account only those cases which require a single 
asymmetry measurement (see Table X of Appendix D). In addition, we consider as observables the cross sections 
and single proton polarizations (xO, yO, zO) only. We have studied these selected structure functions in the same 
kinematic regions as taken for the survey of the unpolarized structure functions but present here only those kinematic 
cases where we found the strongest signatures. 

In Fig. 24 we show the most relevant cases for the IC contribution. It is not surprising that they are all found at an 
excitation energy in the A-resonance region at lower momentum transfer. As pointed out before, appropriate potential 
models with A degrees of freedom are lacking and thus we only show results for the coupled channel calculation with 
the renormalized BONN-Qb potential. The figure shows that there are quite a number of different structure functions 
which are essentially dominated by their IC contribution and thus are ideal cases for the study of such IC effects. It 
it worthwhile to note that there are not only transverse (/^ 10 , f^ 1 , /t°( z 0), /t 1 ^))) but also LT-type structure 
functions (/^"(zO), /^ + (a;0)) among the cases presented. 

Observables which exhibit sensitivity to the MEC contribution are shown in Fig. 25. While for the strongest 
relative MEC effect is found around 9 = 90°, for J'- 1 {xQ) and /y?°(xO) one has the most pronounced effects at 
backward angles. The potential model dependence is not very important and is particularly small for /j?°(xO). We 
would like to mention that the result is very similar to that of fjP for the same kinematics (see Fig. 21). In fact, 
/j. is very sensitive to the MEC contribution and can in principle be determined in a single measurement (scattered 
electron at backward angle and <j> = 45°). 

In Fig. 26 we show those structure functions where there are major effects due to relativistic contributions. A 
very strong effect is seen for , although this structure function shows a rather severe potential model dependence 
as well. For /y 11 (zO) one finds a similar situation, but with somewhat smaller relativistic influence and a less pro- 
nounced variation with the potential model. On the contrary, for f LJ T (yO) one observes essentially no potential model 
dependence and a large relativistic contribution at forward angles. 

In view of the fact that, besides the unpolarized differential cross section, the quantities which one measures 
experimentally are asymmetries, we present in Figs. 27 through 29 the asymmetries corresponding to the structure 
functions of Figs. 24 through 26 in order to see which of them produce sizeable, i.e. easily accessible asymmetries. In 
order to emphasize the interaction and relativistic effects in the numerator, the following asymmetries always refer 
to the unpolarized differential cross section in which all contributions from MEC, IC and RC are included, i.e. So in 
Eq. (83) is So(T). 

With respect to IC effects, one readily notes in Fig. 27 three sizeable asymmetries of the order of one, namely A^ d 
for two different settings (9d = 0° and 90° with 4>d = 90°), requiring electron and deuteron vector polarization, and 
A e (zO), requiring electron polarization and a proton polarimeter. Two other asymmetries, A^(zO) and A^ d (zO), are 
of the order of 0.2. Only AY d (xO) is too small being of the order of 0.004. Of the asymmetries sensitive to MEC in 
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FIG. 28: Influence of MEC on asymmetries for kinematic settings for the separation of the various structure functions in Fig. 25 
calculated with the BONN-Qb potential. Notation as in Fig. 27. 
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FIG. 29: Influence of RC on asymmetries for kinematic settings for the separation of the various structure functions in Fig. 26 
calculated with the BONN-Qb potential. Notation as in Fig. 27. 



Fig. 28 the largest one is A v ed being of the order of 0.7, whereas the other two are smaller, of the order of 0.4. With 
respect to RC, two asymmetries in Fig. 29 are sufficiently large, namely A d and A d (yO) with a magnitude of about 
0.1 to 0.2, but only A d exhibits a large RC effect. 

Finally we would like to briefly discuss irrespective of the problem of separating a specific structure function the 
question which observables show in specific asymmetries of sizeable magnitude either a large influence or none by 
interaction contributions. We have investigated for a large number of observables various asymmetries for a variety 
of kinematic settings. Some results of this search are presented in Figs. 30 and 31. 

Fig. 30 exhibits the tensor asymmetry A d for two different settings of the angles <j> 1 Od and <fid- Without interaction 
effects the asymmetry would be quite small, about 0.1 or less. In the left panel, MEC lead to a strong enhancement 
around 0° and 180° which are partially canceled by relativistic contributions. On the other hand, IC show a large 
and dominant influence around 90° and 270° where the other effects are very small. The right panel is an instructive 
example on how the the asymmetry and the relative size of the various contributions change with a change of the 
out-of-plane angle <f) and the deuteron orientation angles. In the left panel of Fig. 31, the vector asymmetry A^ d 
exhibits at low energy but higher momentum transfer a large MEC effect which is cut down to almost one half by 
relativistic contributions whereas IC is negligible. The same asymmetry at a different kinematics, E np — 100 MeV 
and q = 1 fm~ 2 , in the right panel of Fig. 31 shows a strong influence from MEC and IC of almost equal size and 
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FIG. 30: Tensor asymmetry for E np = 250 MeV, q 2 = 1.5 fm 2 calculated with the BONN-Qb potential. Notation as in 
Fig. 27. 
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FIG. 32: Asymmetries at three different on-quasi-free kinematics, calculated with the BONN-Qb potential. Left panel A e for 
E np = 10 MeV, q 2 = 1 fm~ 2 , middle panel A\ for E np = 30 MeV, q 2 = 3 fm~ 2 , and right panel A T d for E np = 250 MeV, 
q 2 = 25 fm~ 2 . Notation as in Fig. 27. 



interfering constructively, while RC effects are tiny. 

As last examples we show in Fig. 32 three asymmetries for on-quasi-free kinematics, two at low energies and one in 
the A-region which exhibit a nearly interaction independent behaviour. The two asymmetries at quite low energies, A e 
in the left panel and in particular in the middle panel show very little influences from interaction and relativistic 
effects. In general these effects increase with growing energy E np . But even for the example in the right panel at 
considerably high energy such influences are still relatively small. It is not surprising, that one finds such a behaviour 
for kinematics belonging to the quasi-free ridge, because one expects that influences from FSI, MEC and IC will be 
minimal there. Because of their model independence, such cases provide consistency checks for theory and experiment. 

With this we will close the discussion of our results. We will not compare to experimental data on the exclusive 
reaction [61-75] because they have already been compared to our approach and thus nothing new can be said. 



VI. SUMMARY AND CONCLUSION 



In this work we have presented a thorough and detailed survey on polarization observables in electrodisintegration 
of the deuteron. It contains a general review of the basic formal aspects with respect to kinematics, the definition of 
observables and structure functions, their multipole decomposition, the question of completeness and the construction 
of independent sets of observables including analytic solutions of the i-matrix elements in terms of structure functions. 

Furthermore, a detailed account of the dynamic ingredients of the theoretical framework is given with respect to 
the basic NN-ioxce, the associated meson exchange currents, isobar configurations for the consideration of internal 
nucleon dynamics, and leading order relativistic contributions. 

As results we first have presented a general survey on the unpolarized and polarized form factors of the inclusive 
process indicating the various kinematic regions with respect to energy and momentum transfer where sensitivities 
to MEC, IC and RC according to the theoretical framework are to be expected. Furthermore, a comparison with 
experimental data for the near threshold region and for quasi-free kinematics is discussed. For the threshold region 
the influence of interaction effects is quantitatively confirmed. Also the quasi-free peak is well reproduced provided 
the reduction by relativistic contributions is included. 

For the exclusive reaction an extensive survey on the unpolarized structure function for a representative grid of 
energy and momentum transfers is given including a detailed discussion of the various interaction effects. Furthermore, 
for a selected set of polarization structure functions whose determination requires only one asymmetry measurement, 
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we have chosen those kinematic regions in which either IC or MEC or RC play a major role with subsequent discussion 
of the associated asymmetries. 

Finally we hope that this survey will stimulate further experimental and theoretical research allowing one to test 
more thoroughly the underlying basic dynamics. This is of particular importance with respect to the question how 
well do we understand the strong interaction in terms of effective degrees of freedom, i.e. in terms of meson, nucleon 
and isobar degrees of freedom, and where will we need to introduce explicitly the basic quark and gluon degrees of 
freedom of QCD. 
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APPENDIX A: EXPLICIT EXPRESSIONS FOR THE U$ XIM IN THE STANDARD REPRESENTATION 

OF THE i-MATRICES 

In view of the angular momentum algebra, it is useful to switch to a spherical representation replacing the cartesian 
components of the nucleon spin operators by their spherical ones 

T 

= £ £ 8T a"ri ] ® for « = 0,..., 3, (Af) 

T = 0,1 I/=-T 

where we have introduced the equivalent complete set of 2 x 2-matrices in s = 1/2 space 

a [ J ] (i) forr = 0, 1 and \v\ < r, , (A2) 

defining ctq ' = er = I2 an d a ^ — °V> the spherical components of the Pauli spin matrices. The transformation 
matrix is given by 



s a v = c(a) 5 T .T( a ) (£i/,5?(cO + c(a) S v ,-p( a )) , (A3) 



with 



c(a) = 5 a2 - S a i , c(a) 



1 for a = 0, 3 
for a = 1, 2 



V2 



for a = 0, 3 

1 for a = 1, 2 



(A4) 



r(a) =1-S a0 , v{a) 

For the inverse transformation from spherical to cartesian components one easily finds 

3 

4 T] (<) = £<&, (A5) 

where we have introduced 

C = c{v)(6 aMT . y) + iv 6 aib ( TiV )) , (A6) 

and 

v 

vl' 

Then the transformation of the Ws to spherical components is given by 



c(v) = j= 5\ v \i + 5 v0 , a(r, v) = 3r- 2\v\ , b (r, v) = 3r - \v\ . (A7) 



^ X J M - E s a Ys^U^, (A8) 
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where U*,„,£; is defined as in (51) with a Xj (i) being replaced by the spherical components according to (A2). In the 
coupled representation one has the following explicit form for the W's 



where the quantities Vy^f M are gi ven by 



(A9) 
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Specifying the observable X, one has in detail: 
(i) differential cross section (X = 1 = (00)) 
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(ii) single nucleon polarization (X = (iO) for the proton or X = (Oi) for the neutron) 
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APPENDIX B: QUADRATIC RELATIONS BETWEEN OBSERVABLES 



In this appendix we will show that for a set of n independent t- matrix elements {tj ; j = 1 . . . n} one finds exactly 
(n — l) 2 quadratic relations between observables by which the n 2 linearly independent observables are reduced to a 
set of 2n — 1 independent ones. To this end we introduce the bilinear form in the i-matrix elements 

Tj>j = t*>tj ■ (Bl) 
which can be expressed as a linear form of the observables O a 

If, Xv'/ C '"- ( B2 ) 

a 

with appropriate coefficients r" . They have the property 



- - ry, . (B3) 



a * 

'J'3 jj' ' 



which follows from Tjij = T*y and the fact that the observables are real quantities. It is straightforward to show that 
these bilinear forms obey the relation 

Tj'jTi m = Tj' m Tij , (B4) 

which, expressed in terms of observables, yields quadratic relations between the latter. In particular, choosing k = 
I = m, one finds 

Tj'j = , (B5) 

J-kk 

where k can be chosen arbitrarily. It is also clear that from (B5) one can recover the relation (B4). Thus we only 
need to consider the latter relation, and the question is, how many independent quadratic relations one can find. 

We first note, it is sufficient to consider only one specific k, because from (B5) one can derive straightforwardly the 
analogous relation for any other kl . Second, it is sufficient to consider only the cases f < j, because Tjji = T*,j. The 
remaining relations certainly are independent because of the independency of the i-matrix elements. Choosing then 
first j' = j, the case j' = k yields the identity, whereas for j' ^ k one finds 

Tj'j'Tkk = \Tjik\ 2 , (B6) 

which constitute (n — 1) real quadratic relations 

£ 7ft, r« o a o a ' = T n T n ° a ° a ' ■ (B7) 

aa' aa' 

As next we consider the case i < j for which one has N = n(n — l)/2 different pairs. Again one can discard the cases 
f = k or j = k, because they do not result in quadratic relations, thus ruling out n — 1 relations. Therefore, one finds 
in this case (j' < j) a total number of 

N-{n-l)= l -{n-l)(n-2) (B8) 
different complex quadratic relations of the form 

J2 Tfi Tki O a O a ' = J2 r?k rfi O a O a ' . (B9) 

aa' aa' 

Separating these into real and imaginary parts, one finds as total number of independent real quadratic relations 
between observables 

(n- 1) + 2^(n- l)(n- 2) = (n - l) 2 , (BIO) 

which is just the required number of relations in order to reduce the number of n 2 linearly independent observables 
to n 2 — (n — l) 2 = 2n — 1 independent ones. 
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APPENDIX C: MULTIPOLE EXPANSION OF STRUCTURE FUNCTIONS AND FORM FACTORS 



Here we will list more explicit expressions for the multipole expansion of the structure functions of the differential 
cross section 
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Note, that the tilde indicates the incorporation of the hadronic phase factor e lS ^ in the multipole matrix element. 
Now we specialize further in order to obtain the angular coefficients for the unpolarized differential cross section, 

So = c(fcl ab , 4 ab ) ]T (Wf + prtf] d&iO) + Plt/lt <*-io W cos + pttJtt d K 2 M cos 2^) , (C8) 
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f (/) 00+, K 



by setting I = M = in (C2) through (C7). Writing for simplicity /£y T and f y L ' T ^ TT instead of f L j T and /^^ TT ' J 
respectively, one obtains 



tf = -4*X* £ (-) 

L'n'j'Lixj 



_^'+L +i L> L K\\j> j K\ 
I L £' 1 I 



+ l) ((-f +^ + l^in'j'fjLj) Ke(z s "C L '*(p'f) , (C9) 

P&Vi'W) He(* 5 " WfVj') #i L («)) , (CIO) 

+ l)^°(/iYw)»e(i*"C' I ''*(/iy) JVf (w)) , (Cll) 

P&Vi'W) Re(i 5 » JVf ;*(/xY) #i L («)) • (C12) 



L'n'j'Lfij 



At the end of this appendix, we will give the explicit multipolc decomposition of the various inclusive form factors 
of d(e, e')np, which can be obtained from the (K = 0)-coefficicnts of (C2)-(C7) according to 

F W-M = { _y + M (1 + 6mq) I (/ (,)/M+, _ ^IM-, 0) (C13) 



The unpolarized form factors are given by 



^ = ^E^ i (l^(w)l 2 + I^ i (w)l 2 ). ( C15 ) 

the vector polarization form factors by 

Re^Wwft/y)], (C16) 









{ l ; 


L 






1 


-1 




1 





= 16^ ^ (-)' [ ~ \ n \\" 7 ; K Re[i\tf faff frj)] , (C17) 

LL'jjj, 



and finally the tensor polarization form factors by 



F 20 



= - 16 - 2 /f E f o o o) I l' 1 ^"^^[^'(W)*^^)], (C18) 

S fo -i f) 1 1' i 2 }^ 2Pi ^[^'(^r^(^)], (ci9) 



F 2 ^ 1 = 32tt 2 



- on 

'10 

y 



F |o = 167r2 ^| £ ( _y fa L jA |L' L 2 } e -^te[J\TfW^(w)]> (C20) 



LL'jn 



F 2 TT 2 = le^yl ^ (-F ^ 1 j}c-^»e^f;(wr^M. (C21) 
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These expressions have already been reported before in [1]. The following two additional form factors vanish below 
pion threshold if time reversal invariance holds. However, if one considers the partial inclusive reaction d(e, e')np, 
they become nonvanishing above pion threshold in a coupled channel approach including isobars without explicit 
treatment of the iV7V7r-channcls [1]. 

= -32. 2 V2 J2 (-y ( o \ j) W \ ^-^Sm^'feT^to)], (C22) 

F'lt 1 = -32- 2 /f E ( o ^ tj ^ 2 W^P L W*f(w)]- (C23) 

We would like to emphasize, that if one considers the completely inclusive process d(e,e')X, the corresponding 
additional form factors will vanish as long as time reversal invariance holds. 

APPENDIX D: SEPARATION OF A SPECIFIC STRUCTURE FUNCTION 

The first step in the determination of any structure function is the measurement of that asymmetry to which 
it contributes. This requires already a number of measurements with different electron and deuteron polarization 
parameters, i.e., two for A (X), A e (X), A%(X), and Aj(X) and four for A^ d (X) and Aj d (X). Next one must 
determine how many different settings of the angles 4>, 4>d and 9 d are necessary for the final separation. As one will 
see in the following, the structure functions can be divided into different classes according to the minimum number 
of asymmetry measurements required for their extraction. A quick glance at (83) shows that the six asymmetries 
contain differing numbers of structure functions, taking for example an observable of type A, one finds four in Aq{X), 
one in A e (X), 8 in A d (X), 16 in A d (X), 5 in A^ d (X), and 7 in A^ d (X). Similar numbers are found for an observable 
of type B. 

First we will consider an observable of type A. Obviously, the simplest case is the electron asymmetry A e (X) 
containing only f' LT (X). This means that f' LT (X) can be determined from just one asymmetry measurement. Similarly 
a close inspection of the dependence of the asymmetries on the angles </>, (j> d and 6 d shows that four other structure 
functions and two combinations need only one asymmetry measurement. These are, introducing for convenience 



S V d jJ d {X- 0, 0, 9 d ) = S (<P)A^ d (X; 0, 4>, 8 d ) (Dl) 
in order to exhibit the angular dependence explicitly, 

c(fcl ab , kf>)p LT f™{X) = SX(X;^,t,0), (D2) 

c(k\*\kl? h )p LT f%-(X) = -V25r(X;|,0,|), (D3) 

c(fcl ab , k^ h )p' T f^°(X) = -S£(X;|,0,0), (D4) 

c{k^,k^)p' T f' T l \X) = V2SY d X;(^0,^), (D5) 

c(fci ab , k't) p' LT f£' + (X) = V2SY d (X; \, \, \) , (D6) 

c{k\* h ,kf>)p> LT f™{X) = S e T d (X;|,0,O). (D7) 



Nine other structure functions and seven combinations require only two measurements for their separation. These 
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are 



c(fcl ab , k l 2 * h ) PTT f™(X) 
c(fcl ab , k't) PLT fl^ + {X) 

c(k[* h ,k l 2 * h )p TT f^-(X) 
c(k^ h ,k l i h )p LT ft° T (X) 

c(fci ab , k^)p LT f L }+(X) 
c(k[ ah ,k l ? h )p LT f%-(X) 

c(fcl ab , k l * h )p LT f L 2 T + (X) 
c(fcl ab ,^ ab ) PLT /f T -(X) 

c{k^\k^)p' LT f' L ^-{X) 
c(k\^,k^)p' T f T 2 \X) 
c(k\ ah , k^)p' T f T ^{X) 



c(k\*\k^)p> LT f'£'-(X) 



c(fc> ab , kf>)p' LT f£' + {X) 
c(fc> ab , k^)p' LT f L ^ + (X) 

c(k^,k l i h )p' LT f L ^-(X) 



A^(X;|,0,O)-^(X;|,0,O), 

oV (Y . K V It nVfy. 3?r n 7r - 



2 
3 
1 

V2 
^3 



2^(X;J,0,|)-^(X;| ) 0,| 
SJ(X;^,4>,0)-SJ(X;^,4>,0) 



2 

2 
^3 



Sj(X;n^,e o d )-Sj(X;0,^e° d ) 



Sj(X;*,^)-Sj(*;0,£,«5; 



S£(X;O,^O)-S£(X;|,0,O) 



V2 

_V3 
2 
2 

73 
2 

V3 
2 
3 



SY d (X;0,0,-)-SY d (X;-,0,-) 



Sj d (X;0,^,e d ) + Sj d (X;n,^,e d ) 

^(X;0,|,^)-5j d (X ;7 r,|,^)' 
Sj d (X;|,0,^)-Sj d (X;|,^ 



Sj d (X;-,O,0°) + Sj d (X;-,7r,0°) 
2Sj d (X;|,0,|) + S e T d (X;|,0,0)~ 



2 
2 



(D8) 
(D9) 

(D10) 
(Dll) 

(D12) 

(D13) 

(D14) 

(D15) 
(D16) 
(D17) 

(D18) 

(D19) 

(D20) 

(D21) 

(D22) 
(D23) 



where 6° d = arcos(l/\/3) (see Sect. III). 

The terms Pl/^+Pt/t 1 m A d and PlJi^+Pt/t m can also be determined from two asymmetry measurements, 
i.e., 



C (fcl ab ,^ ab ) (p.fl 1 + PT f¥{x)) = -^[s v Ax-^l) + s v d {x-^,y-) 



c(k\^,k^)(p L ff + p T ff(X)) 



1 



Sj(X;^,t,0) + Sj(X;^,4>,0) 



(D24) 
(D25) 



In order to separate the longitudinal from the transverse part one needs in addition a Rosenbluth analysis. 

Increasing the number of asymmetry measurements to three allows one to determine only one further structure 
function and one combination, namely 



c(fcf b , k^)p TT f^ + (X) 



1 

71 



Sdi x 'i ^ ' 2 ' 2^ + Sd{ x 'i 2' '2") ~ ZSXiX; -, - , -) 



c(fc> ab , k\* h )p TT .f™(X) = 1[SJ(X; l,i,0) + £%(X-,^,4,0)] -SJ(X;^,4>,0). 



(D26) 
(D27) 
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All thirteen vector structure functions are then determined and fifteen of the tensor ones or respective com- 
binations. For the remaining combinations f^ + (X), f^,~(X), f^ + (X), f^,~(X), (pl/i} + PtIt~(X)), and 
(plJi 2 + Prf>P{X)) one needs four measurements in order to determine them. The last two combinations require a 
Rosenbluth separation in addition. 

Now we turn to the separation of structure functions for B-type. Without target and electron polarization one has 
two structure functions, one of which can be determined by one out-of-plane measurements, namely (fj^T(X), <f) = ^), 
and the other needs two, [f^f{X),(j) = f , §)• Electron polarization alone without deuteron polarization leads to 
two other structure functions: f^?°~(X) which can be obtained from one out-of-plane measurement at <j> = -|, and 
fLr(X), which requires two settings. 

For the structure functions with target polarization we find that six of them can be determined by a single asymmetry 
measurement. They are 



c(k\* h ,k l 2 * h )p LT fl° T -(X) 



C(fcl ab , k£ b ) P LTfVf(X) 



C (k 



j ab , k^)p' T f^{X) 
l ab , k't)p' T f T w {X) 




(D28) 



(D29) 



(D30) 



c{k^,k^)p' LT f^ + {X) 
c{k^,k^)p' LT f£+{X) 



c(k- 



(D31) 



(D32) 



-V2S£(jr;|,0,|). 



(D33) 
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Sixteen other structure functions require only two measurements for their separation. These are 

c{k^,k X t)pLTf^ T + {X) 



1 

71 



SX(X^,4>,0)-SJ(X;^,4>,0) 



c{k^,k^) PLT fl}+(X) 
c{k^,k)? h )p LT ff+{X) 
c{k x ? h ,kf h )p LT ff T -{X) 
c(k\ ah ,k l i h )p LT fl 2 +(X) 
c(k\* h ,kl* h )p LT fl 2 T -(X) 

C (fc> ab , kl ah )pTTf^ T -(X) 

c(fc> ab , k't) PTT ff T (X) 

c(k^,k^)p' T f^(X) 

c(k\ ah ,k l i h )p' T f^(X) 
c(fc> ab , kf°)p> LT f'£-{X) 
c{k^,kf°)p> LT f™-{X) 

c(k l ? h ,k^)p' LT f£ + (X) 



s v Ax-,\^\)-s v ax^\] 



S J(X] 7r, -, -) - sJ(x- : o, -, 



5J(X;|,0^)-5J(X;^0,^) 



S^X-Q,y-)-Sl{X-^y-) 

0,|)^Sj(X;|,*,0) 
-SX(X;iy-) + V2SX(X ; y-^). 
SJ(X;^,4>,0) - ±SJ(X;^,4>,0) , 



3 
3 

= -V2 



Sj d (X;0,^,e° d ) + Sj d (X;^,e o d ) 

Sj d (X;0,^,e° d ) + Sj d (X;n,^,e° d ) 
S%d(X; 0, — , — ) - SYd{X; — , — , — ) 



Sj d (X;-,0,O)-S e T d (X; 0,0,0) 



nT/y. ?T 7T 7T T 7T 7T 3?T 

°e<il^ ' 2' 2"' 4^ e <^ ^'^'^T 



c(^ b ,^ b )^ T /i? + (X) = y=; 



c{k^\k't)p' LT f£-{X) 



Sj d (X;0^,e° d )-Sj d (X;n,^,e° d ) 

2sux-,\^)+sux-,\rm 

Sj d (X;0,^,9 d )-S^ d (X;^,9 o d ) 



(D34) 

(D35) 

(D36) 

(D37) 

(D38) 

(D39) 
(D40) 
(D41) 

(D42) 

(D43) 
(D44) 
(D45) 

(D46) 

(D47) 

(D48) 

(D49) 



In this last set of equations as well as for what follows we have made a particular choice of angles to allow the 
separations. This choice is not unique and other choices will lead to different coefficients in the linear combinations. 
The terms pLfj^(X) + ptJt°(X) and pLfj}(X) + pTf^(X) in A d (X) can also be determined from two asymmetry 
measurements, i.e., 



c(fci ab ,^ ab ) PL.rr (x)+ PT f™(x) 



1 



SX(X;^,4>,0) + SX(X;^,4>,0) 



c(k\ ah ,k l i h )p L f 1 L 1 (X)+p T f^(X) = --L[^(X;^0,|) + ^(X;^,0,|) 



(D50) 
(D51) 



In order to separate the longitudinal from the transverse part one needs in addition a Rosenbluth analysis. 

Increasing the number of asymmetry measurements to three allows one to determine only two further structure 
functions, namely 



C (k i r, k^) PTT f^{x) 

c(fci ab ,^ ab )p TT /i° T + (X) 



1 
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SX(X; 1 0, \) + S v d {X- ^,0, \) - 2S v d {X- \, 0, \) 



7T »7T ~\ IT 

SX(X; -, 0, 0) + Sj(X; -, 0, 0) - S v d {X; -, 0, 0) . 



(D52) 
(D53) 
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TABLE X: Number of asymmetry measurements for a structure function of an observable X. The symbol (R) indicates the 
need of an additional Rosenbluth L-T-separation. 



x e A 


IM 


L, T 


LT+ 


LT- 


TT+ 


TT- 


T" 


LT'+ 


LT'- 


10 






1 




2 


1 




2 


11 


2(R) 


2 


1 


3 


2 


1 


1 


2 


20 


2(R) 


2 




3 






1 




21 


4(R) 


2 


2 


4 


4 


2 


2 


2 


22 


4(R) 


2 


2 


4 


4 


2 


2 


2 


X £B 


IM 


L, T 


LT+ 


LT- 


TT+ 


TT- 


T' 


LT'+ 


LT'- 


10 


2(R) 


2 




3 






1 




11 


2(R) 


2 


1 


3 


2 


1 


1 


2 


20 






1 




2 


1 




2 


21 


4(R) 


2 


2 


4 


4 


2 


2 


2 


22 


3(R) 


2 


2 


4 


4 


2 


2 


2 



Three settings also allow a determination of the combination 



c{kf\kT)PLf 2 L \X) + p T f^{X) = 



(*4'4'2 ) + S < ( * ; 4 



4' 2 



, + -^'(X ; |,0,O) 



(D54) 



All fourteen vector structure functions and sixteen of the tensor ones are then determined. For the remaining functions 
/ TT f (A), /y^T(A), f^(X), f?ftf(X) and {pLf"i}{X)+pTf^{X)) one needs four measurements in order to determine 
them. The last combination and the one in (D54) requires an additional Rosenbluth separation. A survey on the 
number of asymmetry measurements needed for the separation of a given structure function is given in Tab. X. 



APPENDIX E: EXPLICIT EXPRESSIONS FOR THE MATRIX REPRESENTATION OF Ux XIM 
For the matrix representation of the U\ XIM in (101) we switch to the spherical representation according to (A8) 

C]^'\X)= ]T s^s™C^«(t'v'tvIM). (El) 

T ! V 1 TV 

The explicit forms of the C™*™^^ V 'tvIM) for the hclicity, hybrid and standard bases are: 
(i) Helicity basis with labeling (X%\ a£\ A<'\ \<p) 

cSlS(^^)=2V3(-)^^-^^f'fl(; / \ 'Vj \ A(* \ A. (E2) 



(ii) Hybrid basis with the labeling (Xp\ X„\ X^\ X^ ) , where for proton, neutron and deuteron the spin projections 
refer to the transverse y-axis, 

c ffiV^/M) = 2V3 Wl 'fc, • (E3) 

where 

= (-y- m 'JT,( 3 , L j )i M '- M d J M > M £). (E4) 

m> \ M ml 2 
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(iii) Standard basis with index labeling (sW, mi' \ \('\ m^) 



C S sZ'x}{t'v'tvIM) = 2V3(-) 1 - Xd+s '-< +T ' +T ss'T , T I I ) 



111 



So- 



2 [t 1 t S \ I s' s S 1 
v' v —a J \ m' —to. — er , 



2 2 ' 

1 i r 

2 2 ' 



} ■ (E5) 



APPENDIX F: EXPLICIT EXPRESSIONS OF THE ONE-BODY CURRENT MATRIX IN PLANE WAVE 
BORN APPROXIMATION IN THE STANDARD REPRESENTATION 

In the plane wave Born approximation (PWBA), the final np-scattering state is replaced by a pure plane wave. 
The deuteron wave function has the form 

(r |lA d ) = (r9cf>\l\ d ) = £ ^l(^|(/l)lA d ) , (Fl) 

and the final state plane wave 



r 

=0,2 



(r \ksm s ) = ( 27r )3 /2 \ sm s) ■ (F2) 

The quantization axis is chosen along k, the relative np momentum in the final state. 

Taking the current in the Dirac-Pauli form, one obtains for the one-body current matrix element in the standard 
representation for the nonrelativistic charge density [32] 

(sm s \p NR (q,P)\l\ d ) = S sl T, C lm*m. F V$j\ lm l)> (F3) 

j = l,2 Imi 

where Ci m d lm denotes a Clebsch-Gordan coefficient (Zm;lm s |lAa)- The leading order relativistic contributions 
arise from the p/M-expansion of the Dirac current and from the wave function boost. The corresponding charge 
density operators are denoted by p R and p B , respectively, 



(sm s \p R (q,P)\l\ d ) = — J] E ^fiK( 2G ^- F y){^^9 2 + ^f}^l^)' ( F4 ) 

j — 1,2 lmim' s 

(sm s \p B (q,P)\lX d ) = ^ E E ^lK^{(-« 2 ^^<+5~' S + 5 27x^)<^l^) 

J=l,2 lmim' s 

-(sjq 2 - U g )S sl S ma m' a (kj\iq- f |Zm ; )| . (F5) 
Similarly, one has for the nonrelativistic transverse current (A = ±1) 

( S m 5 |jf*(g,P)|lA d ) = — E E CKl^^^Kni^+^^r^fcilim,). (F6) 

.3=1,2 lmim' s 

For the leading order relativistic current contributions we again distinguish J R from expanding the Dirac current 
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and two contributions from the boost, J^ x and Jf x , with respect to the nonrelativistic convection and spin current, 



s.A 



pi^4> ( r / 

(sm s \J?(q,P)\l\ d ) = -jg^s E E C ^{M (Fi ^ +Gm ^ 



(sm s \J* x (q,P)\l\ d ) 



(sm s \J^(q,P)\lX d ) 



D iX4> 



j— 1,2 lmim' s 

+c Mj \{k] + <?)s™;< + ^ q s-;;: : <]}{k^i mi ) , 

3=1,2 lmim' s 

-(sjq 2 - H q )5 sl 5 m3m > 3 (kj\iq- f|/m;)| , 



(F7) 



(F8) 



32M 3 



3=1,2 lmim' a 



+ 



((q 2 Sj IL q )S°2f° + fSfif') + 1 2 ^3,xS s iS mam ' s + 2Xkq 2 £ a? (k^)} , (F9) 



j'^3 



where kj = k — sjq/2 with Sj = ±1 for j = 1 and 2, respectively. The FOURIER transform of the orbital components 
of the deuteron wave function is given by 



(p\lmi) = J d 3 re- ii; - F (f\lmi) 

= {-) l / 2 4irYi mi (p) / drrui(r)ji(pr) , 
Jo 



(F10) 



and 



(p \iq ■ f \lmi) = -q ■ V p {p \lm{) . (Fll) 
For convenience the various spin matrix elements have been denoted by 

S1T 3< = <(^)"».ki,m.- mi |(^)lm;> = (-l)- m 'f 3 1 , \)(s\\a^\\l), (F12) 

= e-»+({±±)8m.\28 j iff j ■ (k x q)\{\\)^m' s ) = ^kq(m s - m' 8 )S?7« ro -< m< _ m , {9) , (F13) 

SZf' = ^ X H{\\)^sW 3 x q) x \{\\)lm' s ) = \qS s ™°<di m3 _ m ,{6) , (F14) 

S"' e = e^((~) S m s \2s^ x k) X \{\\)lm' s ) = 2 Sj k(m s m' s )S^- m 'di im ,_ m ,(e) , (F15) 

S£f' = e-^((~)sm s m x (2s j k-q)) x \( 1 - 1 -)lm' s ) = S%'< - S%f- , (P16) 

= ( _ r ™. £ K _ mi) ( * 1 S M 

< ro .-K'W4m»-m i W^7 ,>i <*ll4 11 |l* ,, >- ( F1? ) 

Furthermore, 

U q (0) = 2kqcos0, (F18) 

n hq {9) = ^n,-? 2 , (Fi9) 

U jtX (6) = 2 Sj kd{ fi (6). (F20) 



Taking the current in the Sachs form, the expressions for the nonrelativistic and boost currents are simply obtained 
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from (F3), (F5), (F6), (F8) and (F9) by substituting by G E j. For the other relativistic contributions one finds 
(sm s \p R (q,P)\lX d ) = -g^E E C^ ilm ,\G Ej 8 sl 8 msm ,q 2 + {2G Mj -G Ej )S^ (F21) 

J=l,2 lmim' s 

(sm s \j^?)\i\ d ) = e e CK{^[^+^cr+n M cr] 

— 1,2 lmim' s 

+TL jtX [(G El (2q 2 + fcf) - G M ,q 2 )8 sl 5 msm , - (G Ej - G Mo )S s ^<] }(k j \lm l ) . (F22) 



[1] W. Leidemann, E.L. Tomusiak, H. Arenhovel, Phys. Rev. C 43, 1022 (1991). 

[2] H. Arenhovel, W. Leidemann, E.L. Tomusiak, Phys. Rev. C 46, 455 (1992). 

[3] H. Arenhovel, W. Leidemann, E.L. Tomusiak, Few-Body Syst. 15, 109 (1993). 

[4] H. Arenhovel, W. Leidemann, E.L. Tomusiak, Phys. Rev. C 52, 1232 (1995). 

[5] H. Arenhovel, W. Leidemann, E.L. Tomusiak, Nucl. Phys. 641, 517 (1998). 

[6] H. Arenhovel, W. Leidemann, E.L. Tomusiak, Few-Body Syst. 28, 147 (2000). 

[7] H. Arenhovel, Few-Body Syst. 4, 55 (1988). 

[8] H. Arenhovel and K.-M. Schmitt, Few-Body Syst. 8, 77 (1990). 

[9] V. Dmitrasinovic and F. Gross, Phys. Rev. C 40, 2479 (1989) 

[10] H. Arenhovel, W. Leidemann, E.L. Tomusiak, Eur. Phys. J. A 14, 491 (2002). 

[11] F.R. Halpern, Special Relativity and Quantum Mechanics. Englewood Cliffs, Prentice-Hall 1968 

[12] M.A. Martin and T.D. Spearman, Elementary Particle Theory. Amsterdam, North-Holland 1970 

[13] D.R. Giebink, Phys. Rev. C 32, 502 (1985) 

[14] T.W. Donnelly and A.S. Raskin, Ann. Phys. (N.Y.) 169, 247 (1986). 

[15] E.M. Rose, Elementary Theory of Angular Momentum, Wiley New York 1957. 

[16] B.A. Robson, The Theory of Polarization Phenomena. Oxford, Clarendon Press 1974 

[17] F. Partovi, Ann. Phys. (N.Y.) 27, 79 (1964). 

[18] H. Arenhovel and M. Sanzone, Few-Body Syst. Suppl. 3, 1 (1991). 

[19] M. Kawaguchi, Phys. Rev. Ill, 1314 (1958). 

[20] A. Cambi and B. Mosconi, Phys. Rev. C 26, 2358 (1982). 

[21] A.S. Raskin and T.W. Donnelly, Ann. Phys. (N.Y.) 191, 78 (1989). 

[22] J.M. Blatt and L.C. Biedenharn, Phys. Rev. 86, 399 (1952). 

[23] W. Fabian and H. Arenhovel, Nucl. Phys. A 314, 253 (1979). 

[24] H. Arenhovel, M. Danos, and H.T. Williams, Nucl. Phys. A 314 

[25] W. Leidemann and H. Arenhovel, Nucl. Phys. A 465, 573 (1987). 

[26] F. Ritz, H. Goller, T. Wilbois, and H. Arenhovel, Phys. Rev. C 55, 2214 (1997). 

[27] R. Machleidt, K. Holinde, and Ch. Elster, Phys. Rep. 149, 1 (1987). 

[28] M.M. Nagels, T.A. Rijken, and J.J. de Swart, Phys. Rev. D 17, 768 (1978). 

[29] M. Lacombe, B. Loiseau, J.M. Richard, R. Vinh Mau, J. Cote, P. Pires, and R. de Tourreil, Phys. Rev. C 21, 861 (1980). 

[30] R.W. Wiringa, R. Smith, and T.L. Ainsworth, Phys. Rev. C 29, 1207 (1984). 

[31] R.W. Wiringa, V.G.J. Stoks, and R. Schiavilla, Phys. Rev. C 51, 38 (1995). 

[32] G. Beck and H. Arenhovel, Few-Body Syst. 13, 165 (1992); 15, 37 (1993) (E). 

[33] S. Galster et at, Nucl. Phys. B 32, 221 (1971). 

[34] S. Blatnik, N. Zovko, Acta Phys. Austriaca 39, 62 (1974). 

[35] G. Hohler et al, Nucl. Phys. B 114, 505 (1976). 

[36] M.F. Gari, W. Kriimpelmann, Phys. Lett. B 274, 159 (1992). 

[37] H.W. Hammer, U.-G. Meifiner, and D. Drechsel, Phys. Lett. B 385, 343 (1996). 

[38] H. Arenhovel, Phys. Lett. B 199, 13 (1987). 

[39] E.L. Tomusiak and H. Arenhovel, Phys. Lett. B 206, 187 (1988). 

[40] H. Arenhovel, W. Leidemann, E.L. Tomusiak, Z. Phys. A 331, 123 (1988); A 334, 363 (1989) [Erratum]. 

[41] M. Schwamb and H. Arenhovel, Nucl. Phys. A 696, 556 (2001). 

[42] D.O. Riska, Phys. Scrip. 31, 471 (1985). 

[43] A. Buchmann, W. Leidemann, and H. Arenhovel, Nucl. Phys. A 443, 726 (1985). 

[44] H. Arenhovel and M. Schwamb, Eur. Phys. J. A 12, 207 (2001). 

[45] H. Arenhovel, Z. Phys. A 302, 25 (1981). 

[46] M. Gari and H. Hebach, Phys. Rep. 72, 1 (1981). 

[47] H.J. Weber and H. Arenhovel, Phys. Rep. 36, 279 (1978). 

[48] A.M. Green and M. Sainio, J. Phys. G 8, 1337 (1982). 

[49] W. Leidemann, diploma thesis, Mainz 1980. 



■56 



[50] J.L. Friar, Ann. Phys. (NY) 104, 180 (1977). 

[51] E. Truhlik and J. Adam, Jr., Nucl. Phys. A 492, 529 (1989). 

[52] H. Goller and H. Arenhovel, Few-Body Syst. 13, 117 (1992). 

[53] K. Tamura, T. Niwa, T. Sato, and H. Ohtsubo, Nucl. Phys. A 536, 597 (1992). 

[54] J. Adam, Jr. and H. Arenhovel, Nucl. Phys. A 614, 289 (1997). 

[55] T. Wilbois, G. Beck, and H. Arenhovel, Few-Body Syst. 15, 39 (1993). 

[56] S. Auffret et al, Phys. Rev. Lett. 55, 1362 (1985). 

[57] R.G. Arnold et al, Phys. Rev. C 42, Rl (1990). 

[58] G.G. Simon et al, Nucl. Phys. A 324, 277 (1979). 

[59] B.P. Quinn et al, Phys. Rev. C 37, 1609 (1988). 

[60] B. Boden et al, Nucl. Phys. A 549, 471 (1992). 

[61] H.J. Bulten et al, Phys. Rev. Lett. 74, 4775 (1995). 

[62] A. Pellegrino et al, Phys. Rev. Lett. 78, 4011 (1997). 

[63] K.I. Blomqvist et al, Phys. Lett. B 424, 33 (1998). 

[64] W.-J. Kasdorp et al, Few-Body Syst. 25, 115 (1998). 

[65] S.M. Dolfini et al, Phys. Rev. C 60, 064622 (1999). 

[66] M. Ostrick et al, Phys. Rev. Lett. 83, 276 (1999). 

[67] D.M. Nikolenko et al, Nucl. Phys. A 684, 525c (2001). 

[68] H. Zhu et al, Phys. Rev. Lett. 87, 081801 (2001). 

[69] Z.-L. Zhou et al, Phys. Rev. Lett. 87, 172301 (2001). 

[70] P. von Neuman-Cosel et al, Phys. Rev. Lett. 88, 202304 (2002). 

[71] P.E. Ulmer et al, Phys. Rev. Lett. 89, 062301 (2002). 

[72] U. Miiller et al, Proc. Int. School on Nuclear Physics, Erice, Italy, Prog. Part. Nucl. Phys. 50, 483 (2003). 
[73] D.M. Nikolenko et al, Phys. Rev. Lett. 90, 072501 (2003). 
[74] R. Madey et al, Eur. Phys. J. A 17, 323 (2003). 
[75] R. Madey et al, Phys. Rev. Lett. 91, 122002 (2003). 



